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1 INTRODUCTION 

The emission-line analysis of ionized gaseous nebulae such as 
planetary nebulae (PNe; ionized ejected envelopes of low- to 
intermediate-mass evolved stars) and H II regions (ionized gas 
clouds marking the birth places of stars, like the Orion Nebula) 
is one of the best established methods for obtaining the abundance 
of elements in our own and in external galaxies. Abundance studies 
of PNe reveal the composition of matter returned to the interstel- 
lar medium (ISM) after being processed by stellar nucleosynthesis, 
while those of H II regions provide a snapshot of the composition 
of the ISM in the latest episode of its evolution. 



Photoionization models of the giant H II region 30 Doradus are built and confronted to avail- 
able UV, optical, IR (ISO) and radio spectra, under black-body or CoStar spectral energy dis- 
tributions for the primary source and various density distributions for the nebular gas. Chem- 
ically homogeneous models show very small rms electron temperature fluctuations and fail 
to reproduce the heavy element optical recombination line (ORL) spectrum of the nebula. 
Dual abundance models incorporating small-scale chemical inhomogeneities in the form of 
hydrogen-deficient inclusions which are in pressure balance with the normal composition am- 
bient gas, provide a better fit to the observed heavy element ORLs and other nebular lines, 
while most spectral features are satisfactorily accounted for. The inclusions, whose mass is 
^ 2 per cent of the total gaseous mass, are 2-3 times cooler and denser than the ambient neb- 
ula. Their O/H abundance ratio is ^ 0.9 dex larger than in the normal composition gas and have 
typical mass fractions of X = 0.687, Y = 0.273 and Z = 0.040. Helium is found to be about 
as deficient as hydrogen in the inclusions, while elements heavier than neon, such as sulfur and 
argon, are quite possibly enhanced in proportions similar to oxygen, as indicated by the most 
satisfactory dual abundance model obtained. This suggests that the posited H-deficient inclu- 
sions may have arisen from partial mixing of matter which was nucleosynthetically processed 
in a supernova event with gas of normal LMC composition. The average gaseous abundances 
of the chemically inhomogeneous models are ^ 0.08 dex higher than those of the homoge- 
neous models, yet they are lower by a similar ^ 0.08 dex than those derived from standard 
empirical methods (fully corrected for inaccuracies in ionization correction factors and differ- 
ences in atomic data) which postulate temperature fluctuations in a chemically homogeneous 
medium. Attention is drawn to a bias in the determination of H II region helium abundances 
in the presence of H-deficient inclusions. It is argued that these results provide evidence for 
incomplete small-scale mixing of the interstellar medium (ISM). The case for the existence 
of abundance inhomogeneities in H II regions is examined in the light of current theoretical 
considerations regarding the process of chemical homogenization in the ISM. 

Key Words: ISM: abundances - H II regions - ISM: evolution - ISM: individual: 30 Dor - 
galaxies: individual: LMC 

has arisen from the fact that, for both PNe and H II regions, the 
abundances of C, N, O, and Ne (relative to hydrogen) based on 
optical recombination lines (ORLs) are found to be consistently 
higher than the abundances of the same elements derived from 
their UV, optical and infrared collisionally excited lines (CELs; 
Barker 1984; Peimbert, Storey & Torres-Peimbert 1993; Rola and 
Stasinska 1994; Liu et al. 1995, 2000, 2001, 2004; Esteban et al. 
2005; Tsamis et al. 2003a, 2003b, 2004). The ORL/CEL ionic 
abundance discrepancy factors (ADFs) for O^^ and C^^ have been 
shown to be strongly correlated with the difference between the 
nebular electron temperatures derived from the nebular to auroral 
[O III] optical-line ratio and the H I Balmer recombination contin- 



Since before the 1990's, a source of controversy in the field 
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uum break at 3646-A (Liu et al. 2001; Tsamis 2002; Tsamis et al. 
2004). 

Deciphering the thermal and density structure of ionized neb- 
ulae is a very important matter in its own right, given the exponen- 
tial dependence of CEL-based abundances on the adopted nebular 
electron temperature {Te, e.g. Osterbrock 1989) and the possibil- 
ity that such abimdances might be largely in error, i.e. imderesti- 
mated, in the case of strong fluctuations in temperature (e.g. Pe- 
imbert 1967) or density (e.g. Rubin 1989) within the nebular vol- 
umes. HST observations of the optical [O III] nebular to auroral 
line ratio in the Orion nebula, where the ADF is small, revealed 
residual small scale temperature fluctuations (O'Dell et al. 2003), 
whose origin has not been fully elucidated. However, it has now 
become clear that invoking in all generality the presence of such 
fluctuations in a chemically homogeneous medium fails to account 
consistently for the observed behaviour of the ADFs; chiefly the 
fact that these do not correlate with the excitation energies or crit- 
ical densities of the lines in question. This is true for PNe (Liu et 
al. 2000; Tsamis et al. 2003b, 2004) and for H II regions (Tsamis et 
al. 2003a, hereafter T03). In parallel, the recent downward revision 
of the solar oxygen abundance (Allende Prieto et al. 2001) has re- 
moved the long-lasting impression that the oxygen abundances of 
galactic PNe were lower than the Sun's and possibly therefore only 
lower limits: mean PN oxygen abundances derived from forbidden 
lines (Kingsburgh and Barlow 1994; Exter et al. 2004) are now in 
good agreement with solar values. Finally, recent spectroscopic re- 
sults making use of the weakly Tc -dependent ratios of O ii ORLs 
have provided evidence for the presence of a very low temperature 
(< 10^ K) gaseous component embedded within both typical PNe 
(Tsamis et al. 2004) and at3fpical ('bom-again') PNe (Wesson et al. 
2003). 

The question of whether the supersonic velocities and in- 
duced shocks from Wolf-Rayet stellar winds and SNRs (e.g. Chu 
& Kennicutt 1994) might produce significant temperature fluctua- 
tions within giant H ii regions (GHIIRs), tmaccounted for by static 
photoionization codes, warrants further analysis: however, as long 
as the ORL/CEL abundance discrepancies are found to be uncor- 
related with the excitation energies of the various collisionally ex- 
cited lines, the answer to the above question does not bear heavily 
on the ORL/CEL abundance discrepancy issue. 

The currently most viable explanation for both the ORL/CEL 
abundance dichotomy and the very low temperatures derived from 
recombination lines or continua, posits the existence within neb- 
ulae of a gas component enhanced in heavy elements and there- 
fore efficiently cooled by emission in fine structure CELs (Liu et 
al. 2000; Pequignot et al. 2002; Tsamis et al. 2004). This mate- 
rial which is traced by the emission of weak heavy element ORLs, 
is intermingled with a more normal composition, hotter (~ 10* K) 
medium, traced by the emission of strong CELs. The above pic- 
ture has gained in credibility since it has been supported by dual- 
abundance photoionization models of PNe (Pequignot et al. 2002, 
2003; Tylenda 2003; Ercolano et al. 2003; Tsamis & Pequignot 
2005 and to be submitted). 

From the standpoint of chemical abundance analyses, an im- 
portant result of the dual-abundance model is that the high metal- 
Ucity component comprises only a small fraction (~ 1 per cent) 
of the total ionized nebular mass in the objects considered so far, 
so that the heavy-element content of the nebulae is closer to the 
one given by the classical forbidden-line analysis (Pequignot et al. 
2002; Tsamis and Pequignot 2004). This contrasts with the clas- 
sical temperature fluctuation paradigm, which invariably leads to 
CNONe abundances closer to those given by the ORL method 



for PNe (Ruiz et al. 2003; Peimbert et al. 2004) and for H II re- 
gions (Peimbert 2003, hereafter P03; Esteban et al. 2005 and refer- 
ences therein). Since emission-line studies of H II regions provide 
a wealth of information regarding elemental abundances in extra- 
galactic systems and therefore have a considerable bearing on cos- 
mic chemical evolution studies, the question of which description 
provides the right answer is a pressing one. 

The co-existence, in PNe, of gas phases with very different 
compositions and temperatures may plausibly result from the ejec- 
tion of material from the evolving central star whose surface com- 
position changes as nuclear reactions and mass loss proceed, as in 
e.g. the 'bom-again PN' scenario (Iben et al. 1983). This is per- 
haps more unexpected in H II regions which are recently ionized 
zones of the 'normal' interstellar medium (ISM). These objects, 
however, are not bom ex nihilo. Their chemical composition de- 
pends on the past history of the host galaxy: the thorough mixing 
of the ISM is a long process in which successive generations of 
massive stars may be involved (e.g. Tenorio-Tagle 1996). They are 
also dynamic systems, hosting potential sources of chemical inho- 
mogeneities, such as supernova remnants (SNRs), evolved massive 
O/WR-type stars subject to extensive mass-loss via stellar winds, 
and photo-evaporating protoplanetary disks. In the face of recent 
studies which have uncovered fairly large ADFs (~2~5 for O'^"'') 
for galactic and extragalactic H II regions (Esteban et al. 2002; 
T03), efforts should be directed at clarifying the situation. 

In this paper we present a photoionization model study of the 
30 Doradus nebula for which high quality, multi-waveband spec- 
troscopy has become available in recent years, including the heavy 
element ORLs of interest. We would like to investigate whether, 
by analogy with PNe, the existence of chemical inhomogeneities 
could provide a solution to the ORL/CEL problem for this class of 
nebulae as well, and thus help us choose between the two conflict- 
ing abundance sets. In Section 2 we collate available spectroscopic 
data pertinent to this study, while in Section 3 we present the mod- 
els and results. Finally, in Section 4 we discuss the implications 
of this work for emission-line studies of photoionized nebulae and 
consider various hypotheses on the likely nature and origin of the 
postulated chemical inhomogeneities in H II regions. 



2 OBSERVATIONS 

At an adopted distance of 50kpc, the 30 Doradus nebula in the 
LMC is the closest extragalactic giant H II region in our cosmic 
neighbourhood. Its full extent in Ha emission is well over a degree 
in diameter (~1 kpc). It is photoionized by a tremendous amount 
of Lyman (Ly) continuum radiation whose total flux has been esti- 
mated to be up to '^lO®^ photons s~^ from radio observations of the 
entire nebula (Mills, Turtle & Watkinson 1978), while almost half 
of this (4.2x 10^^ photons s~^) arises from within the inner lOpc, 
provided by 117 massive OB and WR-tj^je stars of the R136 clus- 
ter in the nebula's core (Crowther & Dessart 1998). The integrated 
(reddened) Ha flux of the central 30' is 1.71x10^* ergcm^'^ s~^ 
(Kennicutt & Hodge 1986), corresponding to 4.2 x 10®^ photons 
s~^ in the Ly continuum. 

Over the years the nebula has been observed extensively 

1 30 Dor contains both SNRs (such as N157B; e.g. Lazendic et al. 2003) 
and numerous Wolf-Rayet stars (see Parker 1993; Crowther & Dessart 
1998), while a plethora of 'proplyds' have been identified in e.g. the Orion 
nebula. 
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from UV to radio wavelengths. In this study we made use of op- 
tical emission-line fluxes by P03 (UVES/VLT 8.2-m) and T03 
(EMMI/NTT 3.6-m); UV fluxes by Dufour et al. (1982; lUE) and 
Gamett et al. (1995; HST FOS); IR fluxes by Vermeij et al. (2002; 
ISO SWS & LWS), except for the [Ar II] 6.98 ^m line noted by 
Giveon et al. (2002); and radio fluxes by Filipovic et al. (1995). 
A compilation of these is provided in Table 3 [col. 3; dereddened 
intensities expressed in units such that 7(H/3)=10*]. 

The ISO SWS and LWS data are for the position '30Dor#3' 
of Vermeij et al. (2002); the positions observed optically by T03 
and P03 fell within this pointing and the 80" LWS circular beam. 
Peimbert's 3" x 10" echelle slit was placed on a bright filamen- 
tary 'ridge' located 64" north and 60" east from HD 38268 (= 
R136) and was contained within the largest SWS aperture, while 
the long-slit of T03 fell about 20-30" to the west, missing the 
20"x33" SWS aperture and passing through the multiple stars 
R139 and R140^ on an approximate NE-SW direction. From a 
comparison of the intensities of 8 H I lines (Bra, Br/? and Pfa-Q 
listed by Vermeij et al. with the theoretical recombination-line ra- 
tios by Storey & Hummer (1995), the estimated H/3 flux, corrected 
for interstellar extinction, falling into the 14"x20" SWS aperture 
is 1.27x10"^° ergcm"^ s"\ The line fluxes between 12-180 /xm 
from Vermeij et al. were scaled down by factors of 1.35 (SWS, 12- 
27.5 ^m), 2.35 (SWS, 27.5^5.2 ^m) and 7.35 (LWS, 45-180 ^m), 
based on the relative sizes of the four rectangular SWS apertures 
and the need for flux continuity in the spectral continuum, and then 
cast in units relative to H/3. Regarding UV data, the dereddened 
intensity of C III] A 1909 is from observations in the years 1979- 
1980 (Dufour et al. 1982), and was adopted instead of more recent 
HST observations by Gamett et al. (1995a) taken through the 1" 
FOS circular aperture. The 10" x20" elUptical aperture of the lUE 
was centered on a position roughly identical to that targeted by the 
SWS. We used the measured line strengths, relative to A1909, by 
Garnett et al. (1995a) for the O III] A1666 and N ill] A1750 inten- 
sities (which are therefore tentative only), while the Si ill] AA1883, 
1892 fluxes were from the F(Si iii])/F(C III]) ratio measured by 
Garnett et al. (1995b) with the FOS. 

The heavy element ORLs utilized in this work include the C II 
A4267 line and the multiplets N ii A5679 (V3), O I A7773 (VI), 
O II A4075 (VIO), O II A4340 (V2) and O ii A4650 (VI) (co- 
added intensities; standard theoretical relative line ratios were used 
to estimate unobserved multiplet components where necessary; see 
e.g. Liu et al. 1995, 2000; Tsamis et al. 2004). Upper limits were 
adopted for the undetected N ii A4041, Ne ii A4392 and Mg ii 
A4481 lines. 

For this study the adopted total dereddened H/3 intensity from 
30 Dor is 1.19x 10~® ergcm"^ s~^, in line with the aforemen- 
tioned Ha estimate by Kennicutt and Hodge (1986) for the central 
30'; the adopted logarithmic reddening constant was c(H/3) = 0.44. 
T03 and Vermeij et al. (2002) found values of 0.41 and 0.48 from 
fixed long-slit and drift-scan slit spectra respectively. These esti- 
mates compare favourably with a mean extinction of about 0.5 dex 
derived from Ha/H/3 ratio maps of the central 6'x6' region by 
Lazendic et al. (2003). The optical line intensities listed in Table 
3 are essentially taken from P03, who dereddened the fluxes us- 
ing c(H/3) = 0.92 (the maps of Lazendic et al. do indicate higher 
than average extinction at the position observed by P03). The NTT 
long-slit spectra of T03 which sampled a more extended 'slice' of 



^ Both R139 and R140 contain WR components; see their spectrograms by 
T03. 



the nebula were consulted for consistency checks. The UV, IR and 
radio fluxes, although carefully calibrated relative to H/3, do not 
necessarily sample exactly the same emission regions targeted by 
the optical observations. 

3 PHOTOIONIZATION MODELLING 

3.1 The code 

The modelling is undertaken using the detailed photoionization 
code NEBU (e.g. Pequignot et al. 2001) in which the radiative trans- 
fer is computed using the outward only approximation along 20 di- 
rections in spherical symmetry, and the ionization state and plasma 
temperature are determined by solving the relevant equilibrium 
equations. A model is built by combining a user-defined number 
of sectors extracted from different spherically symmetric models 
sharing the same central source of primary ionizing radiation, each 
with its own covering factor a;/47r of the source, radial hydrogen 
density distribution N-aij-), gas filling factor e, and set of elemen- 
tal abundances relative to H. The code version used here does not 
include an interstellar dust component. 

3.2 Modelling strategy 

3.2.1 The 'average filament' 

30 Dor is much larger than R136 and can be viewed as a gaseous 
shell photoionized by a point-like source located at the centre of 
an inner cavity. Published images reveal the complex, frothy ap- 
pearance of the nebula, whose overlapping loops, arcs, and shells 
are believed to have originated in interacting stellar-wind outflows 
and expanding SNRs (e.g. see the 3- and 6-cm radio and Ha, H/3 
images by Lazendic, Dickel, & Jones 2003 and references to imag- 
ing and kinematical studies therein). The mean square root electron 
density (A'e^) = of 30 Dor, derived from the H/3 flux of the central 
30' under the assumption of a uniformly filled sphere, is two orders 
of magnitude smaller than the A^c derived from the [O II] and [S ll] 
doublet ratios. This indicates that the bulk of the emission arises 
from a small fraction of the total volume, as suggested also by the 
highly filamentary appearance of the nebula. This common feature 
of GHIIRs is often handled with the concept of a 'filling factor', e. 
In the classical form of this concept, the emitting gas is assumed 
to belong to infinitesimal, optically thin clumps distributed more or 
less uniformly between an inner and an outer radius, and occupy- 
ing altogether a fraction e of the nebular volume (e.g. Osterbrock 
1989). In this approximation, which aims at emphasizing the ob- 
served geometrical extent of the H II region, the high excitation 
gas generally lies close to the energy source whereas the lower ex- 
citation gas is found farther out. This is usually not borne out by 
observation: although shaUow radial gradients may exist through- 
out the nebula, filaments tend to emit both high and low ionization 
lines at whatever distance from the ionizing source.'"* The very ex- 
istence of a filamentary structure contradicts the classical concept. 
As confirmed by photoionization modelling, given the primary ion- 
izing flux and the gas density of 30 Dor, matter becomes optically 
thick to ionizing radiation over distances which are much smaller 
than the nebula and similar to the typical thickness of the observed 

^ This is indicated from an analysis along ~150" of the sUt of T03 which 
was placed roughly along the nebula's radial direction and crossed several 
filamentary nebulosities (cf. fig. 5 of that paper). 
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Table 1. Model parameters and observational constraints. 

Parameter Constraint (flux or flux ratio) 

c(H/3) Ha/Hj9 (taken from P03) 

L* absolute flux 7(H/9) 

Tcff [O in] A4959 + A 5007 

5 (4ryd)'' He II A4686 (for blackbody models only) 

Rin [O III] (A4959 + A5007)/[O II] (A3726 + A3729) 

Pin [OlII]52/im/88/imor[ArIV]/[ArIII] 

Pout [O II] A3726/A3729 

He mean of He I AA4471, 5876 

C C III] A1908 

N [N II] A6548 + A6584 

O [O III] A4363/(A4959 + A5007) 

Ne [Ne iii] A3869 + A3968 

Mg mean of Mg l] AA4562, 457 1 

Si mean of Si ill] AA1883, 1892 

S mean of [S iii] A A63 12, 9069, 9532 

CI mean of [CI iii] A A55 17, 5537 

Ar [Ar III] A7135 + A7751 

Fe [Fe III] A4658 

" Scaling factor of the blackbody flux at >4 ryd (see text). 



filaments (e.g., compare the spatial thickness of the model nebula 
with its radius in Table 2). Large-scale radial gradients, if present, 
are likely due to a global variation of the ionizing parameter with 
position, e.g. in the case that the average density of the different 
filaments does not scale as the inverse square of their distance to 
the ionizing source. 

As an alternative to the classical description, we consider 
macroscopic, essentially radiation-bounded filaments which are 
distributed throughout the nebula. Each one of them produces a 
radially extended shadow, which emits much less than the filament 
itself since it is only subjected to the weak, extremely soft, dif- 
fuse field from other filaments. Here, this alternative description 
is adopted in its simplest first approximation, namely in terms of 
one representative 'average filament'. Given its small geometrical 
thickness, the radial structure of the filament can be treated assum- 
ing spherical symmetry without loss of generality. As a matter of 
fact, the adopted geometry enters only in the computation of the 
diffuse field, and there only to a limited extent since the diffuse 
field is either relatively unimportant in the limiting case of small 
opacity/optical depth, or very local (i.e. fully independent of geom- 
etry) in the opposite limiting case. Depending on the overall cov- 
ering factor of the entire set of filaments, ionizing radiation may or 
may not leak out from the nebula. In this work, it is assumed that 
the overall covering factor is unity, and the power of the ionizing 
source is the one inferred from the integrated Ha emission over 
30', i.e. a large fraction of the total emission. This amounts to com- 
puting the emission of a fully spherical, geometrically thin shell 
with no escape of the primary and diffuse field. Adopting the spec- 
trum secured by P03 as representative of the whole nebula simply 
means that the properties of our 'average filament' will be some- 
how weighted by the regions effectively observed. 



3.2.2 Model parameters 

Computations are therefore done assvmiing a smooth small-scale 
density distribution, i.e. a gas filUng factor of unity. The pair of (Tg, 
A'e) varies according to the following general law for a variable gas 



pressure P, given here as a function of the radial optical depth, t, 
at 13.6 eV, 



n( \ -fin H" -fout , -fout -fin , — 1 

P(t) = -I- tan 



K log( — 



.(1) 



The pressure is related to the gas temperature and density via the 
ideal gas law. At the first step of the computation (r = 0) the ini- 
tial pressure is Pin, while at the last step (t = oo; in practice »1) 
the final pressure is Pout. Eq. 1 introduces 3 free parameters in 
the modelling process: Pn, Pout, and the optical depth at which 
the transition from inner to outer pressures occurs (t„). The 'slope' 
of this transition is controlled by k, which is arbitrarily taken as k 
= 10, ensuring a smooth transition in all computations presented 
here. In practice, will have almost the same value 3.4 ± 0.1 in 
all computations and the gas distribution is practically controlled 
by two free parameters, Pn and Pout, along with the inner neb- 
ular radius. The realistic picture of a filament core surrounded by 
a dilute halo or envelope dictates that Pn < Pout. The adopted 
pressure law proved to be a convenient tool for the exploration of 
the effects of the density distribution on the model constraints and 
predictions; e.g. Pn could be constrained by the [O ill] 52-/im/88- 
/xm line ratio and/or the [O iii]A5007/52-/L(m ratio, or else some 
ionization-balance consideration (e.g. [Ar IV]/[Ar III]), while Pout 
was constrained by the [O ii] A3726/A3729 ratio. 

Either a black body (BB) or a model stellar atmosphere spec- 
trum were used to simulate the spectral energy distribution (SED) 
of the primary ionizing source. The CoStar model stellar atmo- 
spheres (Schaerer & de Koter 1997), which include non-LTE treat- 
ment, line blanketing and stellar wind effects, were adopted. The 
central source was controlled by a luminosity L and a temperature 
Tcft . For the BB models, a scaUng factor 5 was appUed to the flux at 
>4 ryd in order to fit He ii A4686 (from T03). Given the very small 
amount of radiative energy involved here, 5 is of negligible conse- 
quence for any other output or prediction of the model (6 ~ 0.10 
in adopted models). For the CoStar models, additional parameters 
would be (i) the metallicity of the stellar atmosphere, taken to be the 
metallicity of the LMC, and (ii) the luminosity distribution of stars 
in the central cluster. Since our aim was not to realistically model 
any particular stellar cluster but only to evaluate the consequences 
of varying the primary SED within reasonable limits (i.e. BB vs. 
CoStar), stars with luminosity of 10^® ergs~^ were systematically 
adopted here, with their number defining the total luminosity of the 
source. 

In summary, there are 17-18 basic line intensities and inten- 
sity ratios which constrain an equal number of model parameters: 
i.e. the reddening correction; 2 or 3 central source parameters (L, 
Tcft, S); 3 nebular parameters (Pin, Pn, Pout); and 11 elemen- 
tal abundances relative to H - see Table 1. All other UV, optical 
and IR observables (Table 3) were treated as 'model predictions', 
whose relevance was evaluated in terms of the employed atomic 
data, the observational (and calibration) uncertainties, as well as 
more general astrophysical considerations. 



3.3 Chemically homogeneous models 

We initially constructed single-component, chemically homoge- 
neous models. Adopted parameters and derived physical conditions 
for two of these - referred to as SI and S2 (with a BB and a CoStar 
ionizing spectrum respectively) - are listed in Table 2. Chemi- 
cally inhomogeneous, dual-component models Dl and D2 are in- 
troduced in Sect. 3.4. The models are radiation bounded (optical 
depths at the end of each computation are given in Table 2). CoStar 
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models have a larger inner radius and a smaller spatial thickness, 
driven by the need to lower the ionization state of the gas which was 
otherwise overestimated, as judged from the ionization balance of 
oxygen mainly, due to the increased hardness of the primary spec- 
trum. The density profile of model S 1 is quite similar to the profile 
of the normal abundance component of model Dl (Sect. 3.4) plot- 
ted in Fig. 1. The ratios of predicted over observed line intensities 
(hereafter the 'departure ratios') are listed in Table 3 (cols. 4 and 
5) and the overall quality of the fit (i.e. closeness to unity) can be 
judged from those. 

3.3.1 Comments 

Models SI and S2 provide us with a view of the influence of 
the primary spectrum on the predictions. To first order these 
models were successful in that it was possible to generate a set of 
parameters fulfilling all basic requirements of Table 2, while many 
predictions compare well with observation (Table 3). Many of the 
unsatisfactory departure ratios correspond to poorly determined 
line intensities or upper limits. This overall agreement suggests 
that basic assumptions are quite acceptable, while the atomic 
data and the observations used are generally quite accurate. A 
commentary on remaining discrepancies is of interest in order 
to evaluate their significance and determine the astrophysical 
relevance of the chemically homogeneous description. Also, many 
features of Sl-2 are common with the following, more elaborate, 
models. Optical lines are considered first. 

He I. The A7065 line appears overestimated, but can be affected by 
telluric absorption: this line is 1.5 times stronger in the spectrum 
of T03. Singlet He I lines were computed assuming case B recom- 
bination (e.g. Osterbrock 1989). Except for the very weak A4437 
line, all of those singlet lines which are much weaker under case A 
(they are noted with a 'B' in col. 3 of Table 3) are systematically 
overestimated by the model, indicating that there is some departure 
from case B in 30 Dor. This is possible even in optically thick 
nebulae since the resonance lines of He I can be absorbed by 
hydrogen atoms before being converted into optical He I lines, and 
this misfit has no bearing on either the quality of the observations 
or the relevance of the model premises. The He i A4121 misfit 
is possibly due to the observed line intensity being affected by 
blending with O II V20 multiplet hnes other than A4119.22 Usted 
by P03. The perfect fit to A3888 was obtained partly thanks to a 
secondary parameter, namely an average 'He"*" turbulent tempera- 
ture', which is introduced to modulate somewhat the optical depth 
in the 2^S - n'^P transitions, and is proportional to the column 
density of the 2'^S metastable level (concerning the He I triplet 
transfer, see, e.g. Osterbrock 1989). The validity of this description 
and the quality of the observations can be appreciated from (i) 
the good fit to A3 187, also sensitive to optical depth, and (ii) the 
fact that the adopted He+ turbulent temperatures (6, 3, 4, and 
5 X 10* K for SI, S2, Dl, and D2 respectively) are very reasonable. 

Low-ionization lines. Inspection of Table 3 shows that emission 
from low-ionization species, such as [O I] AA6300, 6363, is rather 
well fitted even though model S2 overestimates the [O I] flux by 
~0.3 dex. An overestimation of [O I] may suggest that, contrary to 
our assumption, the material is not strictly radiation-bounded along 
all radial directions. It may also indicate that too many 'hard' pho- 
tons arise from the source and are eventually absorbed at the ioniza- 
tion front (see also the commentary on [Ar iv] below). On the other 
hand, emission of the [N I] A5200 multiplet and the [S ii] lines is 



underestimated by 0.7 and 0.3 dex respectively in model S I . CoStar 
model S2 reproduces the [S II] flux much better, without improving 
the fit to the [N I] lines. Note the overall agreement of the [S II] rel- 
ative line intensities. The intensity of the far-red [S ii] multiplet was 
based on the strongest line [S II] A 10320, re-calibrated by means 
of the 3 strongest neighbouring H I and He I lines, which consis- 
tently point to an upward correction of the intensities by a factor of 
1.28 in this wavelength range of the VLT spectrum (in this range, 
lines with dereddened intensities less than 0.1 in the present units 
lead to inconsistent results). The collision strengths for [N I] are 
questionable as recent results by Tayal (2000) inexplicably differ 
from previous ones by as much as one order of magnitude. Here, 
the values obtained by Berrington & Burke (1981) were used. A 
small wavelength interval including [N I] A5200 was unfortunately 
not observed by T03. 

Finally, the [C I] A9850 line is to a large extent excited by 
recombination (~ 50 per cent on average for the whole nebula). 
The fact that it is underpredicted is to be considered in connection 
with the poor fit to the heavy-element ORLs (see below), even 
though the line can also arise from the photodissociation region 
(PDR), which emits lines like the strong [O I] Gi-jim and is not 
modelled here. 

High-ionization lines. All models, S2 especially, overestimate the 
emission of [Ar iv] and probably [CI IV] also (the latter line is 

weak). The recombination coefficients for these species are poorly 
known, and the adopted approximate empirical values were derived 
from an unpublished model of the PN NGC7027 by Pequignot. 
For example, the empirical coefficient is 8 times the radiative 
coefficient in the case of Ar^"*" recombination. The mean fractional 
concentrations of Ar"*"^ and CI"'"'* are moderate (~8 and 13 per cent 
in model SI; --^11 and 21 per cent in model S2), and very sensitive 
to density as these are the highest ions present in the gas. While 
line fluxes by P03 and T03 are generally in fair agreement, the 
[Ar IV] lines are striking exceptions, with the fluxes of T03 about 
70 per cent larger, suggesting that (i) the concentration of Ar^^ 
indeed varies greatly with position in the nebula and consequently, 
(ii) an overestimation of these line fluxes not exceeding ~0.3 dex 
in a model is probably not physically significant. From trials it 
was found that BB models at constant densities (iVe = 300 cm~*) 
overestimate the [Ar IV] and [CI IV] fluxes by ~50 per cent only, 
whereas CoStar models with a positive radial density gradient 
overestimate this flux by up to 0.60 dex (model S2). The latter 
discrepancy would tend to invalidate the prominent continuum flux 
excess (relative to a BB of the same Teff) shown for hv > 40 eV 
by the CoStar stellar atmospheres, in the way they were used 
here. Also, the extended, low density halo of the high-excitation 
zone of our model filament was tailored to fit the [O iii] FIR 
ratio (models Sl-2). The 80" LWS lobe however sampled a large 
area in which the mean density may be predictably lower than 
the one called for by the [Ar IV] UVES ratio; UVES sampled 
a small bright filamentary region whereas LWS targeted an area 
of lower average surface brightness (see fig. 2 of Vermeij et al. 
2002). These discrepancies are therefore understood in terms of 
conflicting density sampling of the IR and optical density-sensitive 
ratios due to the different spatial coverage of each spectral domain, 
and of the uncertainty regarding the exact shape of the primary 
continuum at high energies (see also model D2 in Section 3.4). 

Intermediate-ionization lines. The [S ill] auroral A6312 and neb- 
ular A9069 + A9532 lines are imperfectly matched. The relative in- 
tensity ratio of the nebular lines agrees with theory within 8 per cent 
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(Mendoza 1983), although these lines are often affected by telluric 
absorption; considering this agreement and their large strength, 
more weight was given to them than to A6312 in the determina- 
tion of S/H. Calculations were performed here using the collision 
strengths by Galavis et al. (1995) for [S III], [Ar III] and [CI IV], 
who devoted particular attention to the fine-structure lines and then 
argued that the [Ar iii] results were not controversial (Galavis et 
al. 1998). The colUsion strengths for [S in] may not be of ultimate 
accuracy however: in a trial calculation using the data published by 
Tayal & Gupta (1999), all 3 optical [S III] lines were almost exactly 
matched using the previously determined S/H abundance ratio. 

The [N II] A5755 line is some 10 per cent high in both S 
models, pointing to a possible weakness of these chemically ho- 
mogeneous models (see Sect. 3.4). Here the presumably accurate 
colUsion strengths by Hudson & Bell (2005) were adopted. 

CNO ORLs. The optical recombination lines of C,N and O are 

systematically underpredicted. This is the raison d'etre of the 
present study (Sects. 3.3.2 & 3.4). 

IR lines. The computed fine-structure ratios for the [Ne in] and 
[Ar III] transitions are not sensitive to model assvunptions in the 
conditions of 30 Dor, and the imperfect match of these ratios 
cannot be taken as an inadequacy of the astrophysical assumptions 
themselves. Whatever the cause, these discrepancies suggest that a 
mismatch of ~ 15 per cent for other fine-structure line ratios might 
not be significant. The discrepancy is larger for [Ar III], but [Ar III] 
21.8- jim is one of the weakest IR lines . It is found from both S 1 and 
S2 models that the predicted [Ne ill] 15.6- and [O ill] 51.8-/xm line 
fluxes support a higher Ne/O ratio than their optical counterparts 
do (if only by 15-18 per cent). In this regard, these homogeneous 
models fail to exactly fit the optical and IR Ne and O lines 
simultaneously. In contrast, the Sl-2 models differ significantly 
for [N III] 57-/um, with the CoStar computation leading to much 
more consistent results. The collision strengths used for [N ill] in 
these computations are from Stafford et al. (1994). Were we to use 
the larger values obtained previously by Blum & Pradhan (1992), 
the difficulty met with [N III] in BB models would be exacerbated. 
Nonetheless, the most obvious failure of these models is the under- 
prediction of [Ne ii]. The IR lines are reconsidered in Section 3.4.2. 



3. 3. 2 Two-sector models 

Chemically homogeneous models comprising two sectors were 
also investigated. The idea (which has been tried in the past) was 
to check whether the different mean densities of the two model 
sectors would induce significant global temperature variations and 
whether this would result in a better fit for the CNO ORLs. A sim- 
ple approach was taken by using the same smooth scale density 
law for both sectors, but imposing different densities (by a factor 
of 2) at the innermost point of each. Vastly different starting values 
were not supported by the observed density line ratio diagnostics. 
The model fitting was as good as for model S 1 with the same abun- 
dances for all elements. The fit to the heavy element ORLs did 
not improve at all however (e.g. C ii A4267 still underestimated by 
0.3 dex). The electron temperature difference between the two com- 
ponents amounted to 120K while the mean density differed by a 
factor ~2. Resulting global temperature fluctuation parameters for 
representative ions were t^(H+) = 0.0038 and t^(0^+) = 0.0026. 
Similar low factors were reported by Kingdon & Ferland (1995) 



who investigated chemically homogeneous models incorporating a 
much broader range of density variations. 



3.3.3 Ionizing spectrum characteristics and model radii 

The Ly continuum ionizing flux of model SI is 7.55x10^^ pho- 
tons"^, comparable to the Lyc flux of 4.2x10^'^ photons"^ de- 
rived by Crowther & Dessart (1998) emanating from within a 10 pc 
radius from the massive core cluster R136, or that of ~1.2x lO^'^ 
photon s^^ obtained by Israel & Koomneef (1979) for the inner 
3'x3' (45pcx45pc) of the nebula. By construction, the total H/3 
flux computed by the model fits the measurement of Kennicutt & 
Hodge (1986) with a circular 30' aperture (i.e. a physical radius 
of '^218pc). The absolute radio fluxes are also well accounted for 
since the effective beam sizes used in the radio observations of Fil- 
ipovic et al. (1995) are effectively contained within the above 30' 
aperture. 

Regarding the shape of the ionizing spectrum, the black-body 
temperature of 49.3 kK corresponding to model SI falls within the 
Teff range of 43-51 kK derived by Martin-Hernandez et al. (2002) 
from the known, observed 30 Dor stellar spectral types using the 
Teff-spectral type calibrations of Vacca et al. (1996) and Martins et 
al. (2002). The Tee obtained using the CoStar stellar atmosphere 
models however is 37-38 kK only, again in excellent agreement 
with the findings of Martin Hernandez et al. (2002) who used those 
models and the cloudy photoionization code to derive the stellar 
properties of 30 Dor from the same ISO observations that were used 
in this work. 

The outer radii of these models span a range of 27-35 pc (Ta- 
ble 2), equivalent to 111-144". This is in very good agreement 
with the actual position of the UVES slit at a projected distance 
of 88" from R136. Furthermore, the derived thickness of our 'av- 
erage filament' is 2-3 pc corresponding to 8-12" at 50kpc. The 
filament of model D2 (Section 3.4) is thinner being roughly 1 pc 
wide (see Fig 1). Inspection of the H/J image of the nebula pre- 
sented by Lazendic et al. (2003), shows that the UVES slit fell on 
a prominent filament with an apparent width of about 8". These 
findings suggest that our approach is successful inasmuch as, being 
based exclusively on observed nebular line intensities, the models 
correctly predict both the global ionizing source properties and the 
approximate spatial characteristics - distance from R136 and extent 
- of the region sampled by spectroscopy. 

The gas phase abundances of these chemically homogeneous 
models are discussed later in Section 3.6 along with our preferred 
30 Dor gaseous abundances. Taking into account the above remarks 
and the commentary of Sect. 3.3.1, we find that the majority of 
lines, both CELs from heavy elements as well as Hi and He I 
ORLs, can be reasonably well fitted, with the crucial exception 
of the CNO ORLs: predictions for these were found to be iden- 
tical across models, i.e. consistently underestimating the observed 
ORL fluxes (Table 3). The behaviour of these lines under widely 
different modelling assumptions such as constant or variable gas 
density, and black body or stellar atmosphere ionizing spectrum, 
necessitates a different approach to the modeUing process. This is 
introduced in the following section in the form of a dual abundance 
model nebula comprising two sectors across which the elemental 
abundances (and therefore other properties such as density and tem- 
perature also) vary. 
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Table 2. Model parameters and physical properties. 



Model 


o 1 




rl-poor 


U 1 

Normal 


Action 

iVlcall 


rl-poor 


DA 

Norinal 


iVlCall 


Primary spectrum 


Planck 


CoStar'^ 




Planck 






CoStar'^ 




Eff. temperature, Teg (kK) 


49.3 


37.6 




51.6 






38.0 




Luminosity, L (10* Lq) 


1.12 


1.98 




1.08 






1.90 






7.80 


10.15 




8.25 






10.50 




Thickness (10^® cm) 


8.39 


5.07 




8.23 






2.87 




-Pin (10-10 CGS) 


4.16 


4.16 




4.00 






6.00 




Pout (10-10 CGS) 


12.20 


12.40 




11.1 






11.4 




Covering factor, oj /Air 


1.00 


1.00 


0.080 


0.920 


LOO*" 


0.086 


0.914 


LOO*" 


Filling factor, e 


1.00 


1.00 


0.105 


1.000 


- 


0.072 


1.000 


- 


Inner Af(H+) (cm-3) 


140 


145 


343 


136 


- 


614 


211 


- 


(Afc> (cm-3) 


201 


204 


440 


192 


198 


696 


271 


279 


(r(H+))(K) 


9962 


9927 


5026 


9942 


9716 


3987 


9895 


9654 


t2(H+) 


0.0045 


0.0064 


0.0696 


0.0054 


0.0175 


0.0544 


0.0072 


0.0224 


t2(02+) 


0.0031 


0.0046 


0.0386 


0.0039 


0.0715 


0.0338 


0.0053 


0.0779 


r (13.6eV) 


90 


92 


29 


90 


- 


9 


34 


- 


H/3 fraction 


1.00 


1.00 


0.081 


0.919 


1.00*" 


0.084 


0.916 


l.OO'' 


Gaseous mass (10^ M©) 


151. 


149. 


3.51 


144. 


147.*" 


1.96 


106. 


108.'' 


He/H 


0.0858 


0.0854 


0.1000 


0.0831 


0.0835 


0.1000 


0.0825 


0.0828 


C/H X 10° 


4.33 


4.60 


66.0 


4.63 


5.98 


72.3 


4.96 


6.09 


N/H xlO^ 


1.74 


1.42 


15.3 


1.74 


1.78 


18.8 


1.46 


1.75 


O/H xlO^ 


22.90 


23.20 


200.0 


24.40 


28.30 


200.0 


25.40 


28.30 


Ne/H xlO^ 


4.33 


4.00 


44.0 


4.52 


5.39 


60.0 


4.38 


5.31 


Mg/H xlO'^ 


52.0 


39.5 


200. 


50.8 


54.7 


280. 


40.0 


44.3 


Si/H X 10^ 


25.2 


28.0 


35. 


28.5 


28.6 


217. 


31.0 


34.0 


S/H X 10^ 


61.0 


68.0 


84. 


68.0 


68.4 


455. 


65.0 


71.4 


Cl/H xlO'^ 


0.59 


0.69 


0.8 


0.66 


0.66 


5.0 


0.72 


0.79 


Ar/H xlO'' 


12.0 


12.6 


16. 


13.1 


13.2 


92. 


13.2 


14.5 


Fe/H xlO'' 


11.7 


11.0 


17. 


13.7 


13.8 


81. 


11.6 


12.7 




0.253 


0.252 


0.274 


0.247 


0.248 


0.273 


0.246 


0.246 


Z 


0.005 


0.005 


0.036 


0.004 


0.005 


0.040 


0.004 


0.006 



" Stellar atmosphere model of 10^^ erg s i luminosity and typical LMC metaUicity. 
' Sum of components instead of mean. 



1000 





Ne (hp) 






Te (hp) 








Ne{ri)^^/^^ 






Te (n) 





0.0 0.2 0.4 0.6 0.8 1.0 

Thickness (/2.87e+18 cm) 

Figure 1. Electron density distribution, N^, and temperatures Te (scaled 
down to 0.02Te and 0.14Te respectively for the normal ('n') and H-poor 
('hp') phases) of the dual abundance model D2. 



3.4 Dual component chemically inhomogeneous models 

Models that comprise two components which are intimately inter- 
mixed are now considered. It is asstmied that clumps or filaments 



of relatively H-poor gas are embedded in the main body of the neb- 
ula, which consists of gas of approximately 'normal' LMC com- 
position. It is also assumed that the clumped phase is coexten- 
sive with the H-rich component. This latter assumption is a sim- 
ple guess suggested by evidence from the spatial analysis of T03, 
which revealed a roughly uniform surface brightness ratio of the 
C II A4267 and O ii ORLs to H/3 across a 160 arcsec swathe of 
the inner 30 Dor region. Technically, the results of 2 appropriately 
defined spherically symmetric computations were combined. In or- 
der for a dual abundance model to be obtained, a small sector of 
the normal composition shell was extracted and substituted by a 
sector of the H-poor shell to the effect that the combined emit- 
ted H/3 flux remained almost unaltered. The gas filling factor in 
the main shell was kept at unity, whereas the small filling factor 
of the H-poor shell was adjusted so that the geometrical thickness 
of both shells were identical. The assumption of intimate mixing 
over a scale much smaller than the shell size led to the comple- 
mentary notion of equal pressures for both components at similar 
optical depths. Given the identical geometrical extent of the two 
sectors and the identical amounts of radiation they absorbed and 
reprocessed, the total radiation field experienced by both were ap- 
proximately similar at any optical depth, and the predicted emission 
from the dual-sector model was very close to the one which would 
arise from a single-sector model with many small H-poor clumps 
uniformly embedded in the normal abtmdance gas throughout the 
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Table 3. Comparison of predictions from models SI, S2, Dl, and D2. The observed, dereddened intensities are in units such that 7(H/3)=10^, 
except the radio continuum fluxes whose units are Jy; the model predictions are given as the ratio of predicted over observed values in each 

case. 



Line 


A(A)« 




SI 


S2 


Dl 


LJz 


H, He recombination Unes, optical 


and radio continua 




H/3 


4861.3 


1000': 


1.00 


1 .00 


1 .00 


1 .00 


Ha 


6562.8 


2850 


1.02 


1.02 


1.02 


1.02 


H7 


4340.5 


483 


0.97 


0.97 


0.97 


0.97 


Hel'* 


3187.7 


27.4 


1.05 


1.06 


1.06 


1.03 


He I (B) 


3613.6 


4.02 


1.28 


1.27 


1.24 


1.22 


He I 


3888.6 


51.7 


1.00 


1.00 


1.00 


1.00 


He 1(B) 


3964.7 


8.61 


1.23 


1.22 


1.19 


1.18 


He I 


4026.2 


21.5 


0.94 


0.94 


0.92 


0.92 


He I 


4120.8 


2.60 


0.73 


0.73 


0.70 


0.69 


He I 


4387.9 


5.50 


0.96 


0.96 


0.95 


0.94 


He 1(B) 


4437.6 


0.565 


0.81 


0.81 


0.78 


0.77 


He I 


4471.5 


43.1 


1.01 


1.01 


1.00 


1.00 


He I 


4713.2 


4.81 


1.12 


1.11 


1.06 


1.05 


xie I 


4921.9 


11.31 


1.01 


1.01 


l.UU 


1.00 


He I (b) 


5015.7 


22.04 


1.22 


1.22 


1.19 


1.17 


He I 


5875.7 


121.2 


1.00 


1.00 


1.00 


1.00 


He I 


6678.2 


33.85 


0.99 


0.99 


1.00 


1.00 


He I 


7065.2 


28.30 


1.76 


1.75 


1.62 


1 .00 


He 1(B) 


7281.3 


4.52 


1.42 


1.41 


1.37 


1.35 


Hen 


4686.0 


0.171 


1.00 


1.54 


1.00 


1.88 


BJ/H/3 


3646 


4.53 


0.97 


0.98 


1.03 


1.06 


cont. 


1.40 GHz 


37.6 


1.05 


1.06 


1.03 


1.02 


cont. 


2.45 GHz 


37.0 


1.04 


1.04 


1.02 


1.00 


cont. 


4.75 GHz 


35.8 


1.02 


1.01 


0.99 


n QS 
u.yo 


cont. 


4.85 GHz 


36.0 


1.01 


1.00 


0.98 


0.97 


cont. 


8.55 GHz 


34.8 


0.98 


0.98 


0.96 


0.94 




Heavy-element recombination lines 






[Cl] 


9850.3+ 


0.276 


0.37 


0.34 


0.84 


0.66 


Cll 


4267.2 


0.919 


0.43 


0.42 


1.00 




W TT 

IN 11 


4038.0 


<0.1 


0.53 


0.39 


0.83 


0.86 


N II 


5679.0+ 


0.186 


0.71 


0.52 


1.00 


1.00 


\j 1 


7773.0+ 


0.124: 


0.58 


0.55 


1.58 


1.85 


U II 


4075.0+ 


3.05 


0.62 


0.63 


0.97 


0.97 


\J II 


4341.0+ 


1.54 


0.69 


0.70 


1.05 


1.05 


M TT 
\J 11 


4651.0+ 


3.70 


0.66 


0.68 


1.01 


1.01 


Ne II 


4394.0 


<0.2 


0.44 


0.44 


0.74 


0.65 


Mg II 


4481.2 


<0.1 


0.50 


0.37 


0.61 


0.55 




CoUisionally excited Unes (UV and optical) 




C II] 


2322+ } 


70: 


0.67 


0.86 


0.65 


0.79 


ttti 
[U IIIJ 


+2321 } 










- 


C III] 


1907+ 


150 


1.00 


1.00 


1.00 


1.00 


C IV 


1549 


12: 


0.31 


0.83 


0.37 


0.69 


rw ii 

LIN IJ 


5197.9 


1.29 


0.28 


0.27 


0.33 


0.16 


[Nl] 


5200.2 


0.693 


0.56 


0.53 


0.64 


0.29 


[Nil] 


5754.6 


1.87 


1.08 


1.10 


0.88 


0.98 


[Nil] 


6583.5+ 


152 


1.00 


1.00 


1.00 


l.UU 


Nlll] 


1744 + 


<15 


0.76 


0.59 


0.61 


0.56 


[0 I] 


6300.3+ 


12.5 


1.05 


2.20 


1.37 


1.04 


[0 II] 


3726.1 


586 


1.00 


1.00 


1.00 


1.00 


[Oil] 


3728.8 


637 


1.00 


1.00 


1.00 


1.00 


[Oil] 


7319.8+ 


37.3 


0.93 


0.96 


0.89 


0.96 


Oili] 


1663 


15: 


2.16 


2.19 


2.18 


2.20 


[O III] 


4363.2 


32.1 


1.00 


1.00 


1.00 


1.00 


[O III] 


5006.8+ 


6790 


1.00 


1.00 


1.00 


1.00 


[Ne III] 


3869.1+ 


444 


1.00 


1.00 


1.00 


1.00 


Mgl] 


4562.6 


0.716 


1.05 


1.06 


1.06 


1.05 


Mg I] 


4571.1 


0.617 


0.94 


0.94 


0.94 


0.95 


Si III] 


1882.7 


20.6 


1.04 


1.04 


1.04 


1.05 


Si III] 


1892.0 


14.4 


0.96 


0.96 


0.96 


0.97 
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Table 3. -continued 



Line 




r , 

obs 


SI 


S2 


Dl 


D2 


[Sll] 


4068.6 


8.71 


0.77 


1.37 


0.82 


1.03 


[Sll] 


4076.3 


2.95 


0.74 


1.31 


0.79 


0.98 


[Sll] 


6716.4 


71.2 


0.68 


1.19 


0.74 


0.90 


[Sii] 


6730.8 


65.6 


0.65 


1.13 


0.70 


0.87 


[S II] 


10320.+ 


6.6: 


0.90 


1.60 


0.96 


1.20 


[S III] 


6312.0 


18.1 


0.79 


0.80 


0.78 


0.75 


[S III] 


9530.6+ 


938 


1.08 


1.09 


1.08 


1.09 


[CI II] 


8578.7+ 


0.76: 


1.04 


1.58 


1.08 


1.64 


[CI III] 


5517.7 


4.49 


1.00 


1.00 


1.00 


0.99 


[CI III] 


5537.7 


3.32 


1.00 


1.00 


1.00 


1.01 


[CI III] 


8433.7+ 


0.250 


0.99 


1.00 


0.99 


0.99 


[CI IV] 


8045.6+ 


0.494 


2.58 


4.79 


3.08 


3.92 


[Ar III] 


5191.8 


0.663 


1.38 


1.38 


1.36 


1.34 


[Ar III] 


7135.8+ 


148 


1.00 


1.00 


1.00 


1.00 


[Ariv] 


4711.4 


1.31 


2.55 


3.97 


2.99 


3.14 


[Ariv] 


4740.2 


0.986 


2.52 


3.92 


2.95 


3.13 


[Fe II] 


5158.8 


0.306 


1.12 


1.78 


1.36 


1.39 


[Fe II] 


7155.2 


0.207 


0.78 


1.21 


0.93 


1.07 


[Fe III] 


4008.3 


0.284 


0.81 


0.82 


0.79 


0.81 


[Fe III] 


4658.2 


4.89 


1.00 


1.00 


1.00 


1.00 


[Fe III] 


4701.5 


1.25 


0.98 


0.99 


0.99 


0.99 


[Fe III] 


4733.9 


0.405 


0.88 


0.88 


0.89 


0.89 


[Fe III] 


4881.1 


1.51 


0.74 


0.74 


0.73 


0.74 


rpf ml 
Lrc iiij 


4924.7 


0.341 


0.66 


0.66 


U.D / 


U.D / 


[re IIIJ 


5270.5 


2.44 


1.01 


1.01 


I.UZ 


I.UZ 


FFe TVl 
Lr-c 1 V J 


6734.4 


0.228 


1.72 


1.52 


1.87 


1.44 




Collisionally excited lines (IR) 






IM TTl 

[IN llj 


121.8 Atm 


<2.7 


0.86 


0.86 


1 20 


1 31 


[N III] 


57.3 fjia 


104 


1.51 


1.10 


1.67 


1.36 


[0 III] 


51.8 ijm 


1330 


1.06 


1.08 


1.40 


1.34 


[O III] 


88.3 Atm 


1250 


1.07 


1.08 


1.29 


1.07 


[Ne II] 


12.8 ^tm 


99.0 


0.39 


0.18 


1.05 


0.79 


[Ne III] 


15.6 fim 


650 


0.91 


0.91 


1.22 


1.41 


[Ne III] 


36.0 /xm 


65.0 


0.81 


0.80 


1.07 


1.21 


[S III] 


18.7 (um 


320 


1.05 


1.05 


1.13 


1.53 


[S III] 


33.5 /^m 


550 


0.82 


0.81 


0.87 


1.04 


[S IV] 


10.5 /^m 


600 


1.62 


1.93 


1.83 


1.59 


[Ar II] 


6.98 


2.4: 


0.94 


1.43 


1.03 


1.83 


[Ar III] 


8.99 /xm 


115 


0.95 


0.95 


1.01 


1.27 


[Ar III] 


21.8 ^iin 


11.5 


0.65 


0.65 


0.69 


0.85 



" '+' indicates that predictions are for the total intensity of the multiplet in each case. 

' ':' indicates uncertain flux and '<' an upper limit flux estimate. 

" Fit to dereddened, total H/3 flux /(H/S) = 1.19x 10"** erg s-^ cm-^. 

'(b)' indicates a prediction for case B recombination of some He I singlet Unes (see text). 



whole shell. The pair of (Te, A'e) values varied along the radial di- 
rection throughout both components according to a P(r) law simi- 
lar to that adopted for the single component models (Fig. 1). 

We consider two representative models: Dl which is built un- 
der a BB SED and has fully radiation-bounded components, and D2 
which is built under a CoStar SED. D2 was tentatively fine-tuned 
to fix the [O I] excess of S2, by considering a large, but finite opti- 
cal depth in the normal component (Sect. 3.3.1). The 'one-filament 
approximation' is obviously coarse and the adopted model filament 
must be seen as a best compromise. Specifically, for Pin < Pout/4, 
the [Ar IV] lines become unacceptably strong, even with the proviso 
of Section 3.3.1, while for Pn > Pout, [O ill] 88-Atm and [S ill] 33- 
/^m become too weak compared to [O iii] 52 //m and [S iii] 18 /um 
respectively (see Section 3.4.2). Accordingly, Pn ~ Pout/2 was 



adopted in model D2. Other models with Pn ~ Pout/3 gave more 
weight to the IR [O ill] ratio. Models Dl-2 should be considered in 
comparison to their respective, one-component counterparts, Sl-2. 

3.4.1 Constraining the metallicities of the model components 

For an increasing metallicity of the clumped component, its elec- 
tron temperature decreases due to more efficient cooling from 
CELs while, because of the equal pressures, the electron density 
proportionally increases. The higher density of the H-poor gas 
favours N'^^ and O^^ recombination which produces and O^. 
H"^ and He° are also enhanced there assisting in the charge ex- 
change process with O^"*" and N^"*" which produces extra singly 
ionized species (Fig. 3) and boosts the overall efficiency of O I 
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relative to O II ORL emission. Tlie observed intensity of tiie very 
weak O I ORL being rattier uncertain, an underestimation 

of the line may not be significant, although a substantial overesti- 
mation is. O I A7773 therefore provides an upper limit to the O/H 
abundance ratio for the posited clumped phase. On the other hand, 
for an O/H ratio significantly less than this upper limit, the H-poor 
phase heats up and the colUsional excitation of [N ii] AA6548, 6584 
increases selectively due to the low excitation energy of the lines. 
These two occurrences can help us to bracket the abundances of the 
posited H-deficient regions (see also Section 3.4.2). 

Regarding model Dl, the Si-Fe abundances of the clumped 
component were kept approximately equal by mass to those of the 
normal component, while those of CNONeMg were enhanced. The 
ratio of oxygen mass fraction (relative to H) of the two compo- 
nents is mo/mo = [(O/Kh/iO/HhMnu/n'ir^, where O/H de- 
notes abundance ratio by number and /in is the mass per gram of 
hydrogen. For model Dl, (0/H)i =2x 10~^ was adopted after trials 
and mQ/rriQ = 7.6 was found. Then the CNONeMg abundances in 
the clumps were adjusted iteratively until the observed intensities 
of the ORLs were matched and without violating the adopted upper 
limits in the case of undetected lines. The resulting ratios for the C, 
N, Ne, and Mg mass fractions are 13.1, 9.5, 9.0, and 4.1 (the lat- 
ter arbitrary) respectively. Thus the H-poor phase of Dl has mass 
fractions of X = 0.690, Y = 0.274, md Z = 0.036. 

For model D2 the abundances of Mg-Fe were adopted to be 
7 times larger in the H-poor component than in the 'normal' one; 
this enrichment factor is close to that found for oxygen. Evidence 
in favour of enhanced Si-Fe in the H-poor component is presented 
in Section 3.4.2. The resulting mass fractions of the H-poor phase 
in model D2 are X = 0.687, Y = Q.llZ, and Z = 0.040. 

Considering a more extreme enrichment of the clumped phase 
in heavy elements, a model where the O/H abundance ratio in 
the metal-rich component was 50 times higher than in the normal 
phase of model Dl (with Mg-Ar abundances equal by mass in both 
components), led to much overestimated O i A7773 emission (by 
0.85 dex), but to a Balmer jump approximately equal to that ob- 
served. A variant of this with Mg-Ar abundances also enhanced 
by 50 times relative to those in the normal phase, overestimated the 
Balmer jump by 23 per cent, and resulted in much overestimated to- 
tal O I A7773 and [Ar ii] 6.98 /xm line fluxes - by prohibitive factors 
of ~10 and ~6 respectively. The fraction of total H/9 flux emitted 
from the clumps was 4.6 per cent, the clump r(0^^) dropped to 
978 K, and the mass of the clumped phase was reduced to 254M0. 
The aforementioned over-predictions for the O I and [Ar II] lines 
did not improve in trials where the computation involved a density 
bounded (as opposed to a radiation bounded) H-poor region. We 
conclude that models where the H-deficient phase has an O abun- 
dance larger than about 10 times (by number) that of the ambient 
30 Dor medium are ruled out. 

3.4.2 Model predictions 

The H-deficient phase of Dl-2 is of lower ionization compared to 
the normal phase: there is ~2 times more H" (relative to H+) and 
-^4 times more He" (relative to He+) in the clump medium than 
in the diffuse medium (Table 5). The clumped phase emits only 
about 8 per cent of the total nebular H/J flux, but it eiruts a more 
signiflcant fraction of [N ii] 6548, 6584-A, [N ill] 57 /um, [O ill] 
52-, 88 ^m, [Ne III] 15.6 /im, and [O II] 3727-A; up to --30per 
cent for the optical and IR nitrogen lines and somewhat less for the 
oxygen IR lines in model Dl. In model D2 about 30-35 per cent 
of flie [Ne III] 15.6-, [S iii] 18.7-, and [Ar iii] 8.99-/ttm fluxes are 



emitted from the cool H-poor component because of its enrichment 
in Ne, S and Ar. The flux portion of the C II A4267 line, the O II 
A4650 VI multiplet, and each of the He I AA4471, 5876 and 6678 
lines emitted from the clumped phase is ~58, 32 and 10 per cent of 
the total respectively for both Dl-2. The relative contributions of 
each nebular component to the individual emission line fluxes are 
shown in Fig. 2 for a representative selection of ORLs, optical/UV 
CELs and IR CELs in model D2. 

Other than the CNO ORLs, inspection of Table 3 shows that 
the predictions for [C I] A9850 and [Ne II] 12.8 /im are also im- 
proved by about 0.3-0.4 and 0.4—0.6 dex respectively in models 
Dl-2 compared to the chemically homogeneous models. The for- 
mer line is efficiently excited largely by radiative recombination^ 
of C+ in the clumped phase (see also Section 3.3.1), while emis- 
sion of the latter is favoured by the lower temperature, as weU as 
the 5 times higher concentration of Ne+ there (Table 5). The major 
flux portion of both these lines (~60-65 per cent of the total) arises 
from the H-deficient inclusions. The He I line spectrum is even bet- 
ter fitted in Dl-2 models (concerning lines sensitive to departure 
from case B and A7065, see Sect. 3.3.1). Notably, due to its par- 
tial emission from a cool medium, the predicted flux of the more 
sensitive to collisional excitation He I A4713 line, is improved by 
6 per cent. Main features of the high excitation lines are as in S 1-2 
(Sect. 3.3.1), with Dl somewhat worse than SI and D2 somewhat 
better than S2. On the other hand, D2 is slightly less satisfactory 
than Dl concerning the fit to the Balmer jump. The interesting case 
of [N II] is considered below. 

The IR CELs merit special attention as they can provide infor- 
mation about the enrichment of Ne, S and Ar in the posited H-poor 
component relative to that obtained for CNO; i.e. whether model 
D2 is to be preferred over Dl. Out of 13 IR CELs listed in Table 3, 
[N II] 122- and [Ar ii] 7 are extremely weak (or undetected) 
and their fluxes should be considered uncertain by at least a fac- 
tor of 2. According to Vermeij et al. (2002a), the most uncertain 
fluxes are those of [S iii] 33- and [Ne iii] 36 /itm, which both suffer 
from the uncertain calibration of the SWS-4 detector (by at least 30 
per cent); [Ne II] 12.8 jMm (25-30 per cent accuracy); and the weak 
[Ar III] 22 /im. For the purpose of determining the IR calibration 
self-consistently by comparison to model results, lines which sig- 
nificantly depend on physical conditions should be discarded, and 
several subsets of IR CELs can be identified (Table 4). In addition 
to [S III] 33 /im and, particularly, [O III] 88 /im, whose critical den- 
sities are low, [S IV] 10.5 /im is also strongly model-dependent in 
that, analogously with [Ar iv] (Section 3.3.1), it is a high-ionization 
line sensitive to both the adopted primary spectrum and the density 
distribution. Thus, an accurate calibration could a priori be based 
on 'set IR5', comprising the 5 remaining lines (i.e. IR5: [N III] 57-, 
[O III] 52-, [Ne III] 15.5-, [S III] 18.6-, and [Ar III] 9 pm). Out of 
these [N ill] is peculiar since it has no optical counterpart and its 
predicted flux is again model-dependent. Accordingly, the most re- 
liable, model-independent set is IR4, i.e. IR5 minus [N III] 57 jim. 
A set including all IR CELs with reliable fluxes is IR7, i.e. IR5 with 
flie addition of [O ill] 88- and [S IV] 10.5 /xm. Finally, IRll is the 
set including the 1 1 IR CELs that were well detected. In Table 4 
the means and 1 a scatters in the departure ratios in models Sl-2 
and Dl-2 are provided for the 4 line sets identified above. 

The small scatter for IR4 indicates that an IR calibration 
from accurately measured lines can be reliably defined (with some 



* The collisional contribution to the line excitation is ~ 35 per cent in the 
H-poor component and ~ 50 per cent on average for the whole nebula. 
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Figure 2. The model over observed flux ratios ('departure ratios') for tlie chemically inhomogeneous model D2 (Tables 2. 3): the relative contributions of each 
nebular component to the individual line fluxes are colour-coded. Three major groups of lines are identified. [C I] A9850 is emitted also from the PDR which 
is not modelled here; O I A7773 is very weakly observed; the deficit of [O III] 88- and [S III] 33-/im is presumably due to the lines being also selectively 
emitted by lower density background gas in the 750 fields of view. Amongst the IR lines shown, [Ne II] 12.8- and [S III] 33-^m are less accurately measured. 
The dashed line shows the suggested re-calibration of the IR fluxes relative to the optical (see Section 3.4.2 for details). 



Table 4. Average departure ratio and Icr scatter for IR lines. 



Line set 


SI 


S2 


Dl 


D2 


IR4'' 


0.99 ± 0.07 


1.00 ±0.08 


1.19 ±0.13 


1.39 ±0.08 


IR5'' 


1.10 ±0.22 


1.02 ±0.08 


1.29 ±0.20 


1.38 ±0.07 


IR?': 


1.17 ±0.24 


1.16 ±0.30 


1.36 ±0.21 


1.37 ±0.12 


IRll'* 


0.99 ± 0.36 


0.96 ± 0.44 


1.20 ±0.28 


1.22 ±0.21 



" IR4: [O III] 52-, [Ne III] 15.5-, [S III] 18.6-, and [Ar III] 9 iim. 

IRS: IR4 -I- [N III] 57 Atm. 
" IR7: IR5 + [O III] 88-, [S IV] 10.5 fim. 

IRll: IR7-I- [Nell] 12.8-, [Ne III] 36-, [S III] 33-, and [Ar III] 21.8 /^m. 



reservations for Dl), provided that these Unes are weakly depen- 
dent on physical conditions. This suggests that the employed colli- 
sional/radiative atomic data and ISO observations are of fair accu- 
racy, and that main features of the present models are not grossly in 
error. While our initial IR calibration is suited to models Sl~2, up- 
ward corrections of the IR line intensities by 19 and 39 per cent are 
indicated from Dl and D2 respectively. This shift is not a weakness 
of the D models since the absolute calibration of the IR relative to 
the optical is rather poor: e.g. note that the H/3 flux, as derived from 
the H I IR lines, is more than twice the one estimated by Dufour et 
al. (1982) falling into their /{/£ 10"x20" lobe (comparison made 
after reduction to an equally sized aperture; the SWS and Dufour 



et al. observed a very similar position). Thus an upward calibration 
of the IR lines relative to the optical serves to reduce the disagree- 
ment between these two F(H/3) estimates. In this respect model D2 
should be preferred. 

Upon adding [N III] (IR5), a increases markedly for BB mod- 
els (SI, Dl), whereas it is stable for CoStar models (S2, D2), show- 
ing that, by adopting the latter primary spectrum, the consistency 
between the optical [N II] CELs (used to determine N/H) and the 
IR [N III] line is preserved. Considering IR7, which comprises all 
accurately measured IR CELs, a is rather stable (and large) for S 1 
and Dl, but is dramatically increased for S2 due to the large de- 
parture ratio of [S IV]. In contrast to S2, model D2 is only mildly 
affected thanks to the conjunction of (i) a smaller departure ratio 
for [S IV] due to the relatively high pressure in the inner zone of 
this model, and (ii) an overall enhancement of the IR line intensi- 
ties (the factor of 1.39 upward shift in calibration) partly due to the 
enhanced abundances of Si-Fe in the H-poor gas (see below). 

The full IRl 1 set leads to a general decrease of the averages, 
with a correlatively larger a, partly attributed to the lesser accu- 
racy of the 4 additional lines. Nonetheless, model D2 again fares 
better with a cr of ~20 per cent, probably not much larger than 
the one expected a priori on the basis of the claimed accuracy of 
the employed line fluxes. The decrease of the averages from IR7 
to IRll is mainly due to [Ar III] 21.8- and [Ne II] 12.8 ^im, with 
the former being by far the weakest line of the set and the latter 
strongly model-dependent, and hence relevant for model appraisal. 
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The particularly bad score of S 1 and S2 when all IR CELs are in- 
cluded is to a large extent a consequence of the weakness of the 
predicted [Ne ii] intensity, while the acceptable score of Dl relates 
to an excellent fit to [Ne ii]. Here the dual abundance models con- 
stitute an improvement over the homogeneous models, with a large 
fraction of [Ne II] arising from the H-poor gas. In this respect, D2 
may appear somewhat less successful than Dl; the fact however 
that [Ne ii] is underestimated in D2 does not necessarily imply a 
failure of the model, as the same occurs in an unpublished nebu 
model of NGC 7027 using the same atomic data. Assuming on this 
basis a 50 per cent improvement to the [Ne ii] 12.8 /nm fit and ig- 
noring [Ar III] 21.8 jim, the 1 cr scatter in the departure ratios for 
the full set of IR CELs is only 14 per cent in D2. 

Concerning elements beyond magnesium and owing to the 
large S/H and Ar/H in the H-poor component of model D2, the 
two most representative lines [S ill] 18.7- and [Ar ill] 9/um are 
enhanced relative to Dl, and are in harmony with the most repre- 
sentative IR CELs of N, O and Ne. This is the origin of the remark- 
ably small scatter found for IR CELs in D2 (Fig. 2), indicating that 
elements beyond magnesium are potentially enhanced along with 
oxygen in the posited H-poor gas. 

Another clue to the superiority of model D2 over Dl is pro- 
vided by N II and [N II]. In model Dl, the temperature of the H- 
poor gas is not low enough to fully quench the [N II] A6584 emis- 
sion, as 0/H caimot be very large, given the constraint imposed by 
the weak O I A7773 ORL. Since N/H there is bound to be suffi- 
ciently large to account for the N II A5679 ORL, [N II] A6584 is 
enhanced, which leads in turn to a decrease of N/H in the 'normal' 
component, with the consequence that [N II] A5755, which was 
somewhat overestimated in homogeneous models, is now underes- 
timated (this line cannot be produced in the H-poor gas, because 
of its high excitation temperature). Lowering N/H then leads to an 
increase of N/H in the H-poor gas in order to accotmt for N ii emis- 
sion, thus potentially causing a kind of divergence of N/H in both 
components, with an increasingly bad fit to the optical [N ii] line 
ratio. Although this divergence does not really occur in Dl, it can 
explain why the predicted [N II] A5755/A6584 ratio in Dl is signif- 
icantly low (Table 3), and why the overall N/H in the same model is 
anomalously small (being almost equal to that in S I, contrary to all 
other abundances; see Table 2). On the other hand, the extra cooling 
provided by notably [S III] and [Ar III] in the H-poor component 
of D2 is enough to almost quench [N II] A6584 emission, thus pre- 
venting the 'N II divergence' and leading to an optimal fit for the 
optical [N II] line ratio. In this sense, the dual abimdance models 
are generally more successful than their homogeneous counterparts 
and, among them, D2 is by far the most successful. 

Finally, the scatter in the IR line departure ratios in D2 can be 
scrutinized from yet another standpoint. Among the 4 lines added 
to IR7 to make up IRll, the very strong [S III] 33 /im stands out. 
Adding this line only to IR7, a increases from 12 to 15 per cent. As 
previously noted, [S III] 33- and [O III] 88 /im, are the only lines 
which can be seriously affected by coUisional de-excitation in the 
low density conditions of 30 Dor. The average for the remaining 6 
lines of IR7 is 1.42, close to the re-calibration factor of 1.39 ob- 
tained from IR4, and the scatter is again 8 per cent only (similar 
to that found for IR4 or IRS; Table 4). Thus after re-calibration, all 
strong IR CELs in D2 are predicted in remarkable agreement with 
observation, with the exception of [O ill] 88- and [S ill] 33 /um, 
which are both underestimated by a factor of ^1.3 (Fig. 2), sug- 
gesting that the fields of view of LWS and SWS-4 sampled regions 
of lower average density than the bright 'optical' filament. In a trial 
calculation with Ri^ twice larger and the gas pressure one third that 



of model D2, the IR [O III] predicted line ratio was inverted, while 
the [S III] ratio was better fitted without being inverted; predictions 
for other IR CELs remained fair. This suggests that taking into ac- 
count the background emission falling into the ISO lobe could help 
in improving the fit to [O III] 88- and [S III] 33 jim, although the 
[S III] line ratio would still be off somewhat. At this point we re- 
call that the [S iii] 33 /um line was not amongst the most accurately 
measured, and that the collision strengths for the [S ill] IR CELs 
are probably not of ultimate accuracy; the values by Galavis et al. 
(1995) used in the present investigation, differ very substantially 
from those obtained by Tayall & Gupta (1999). 

In conclusion, the few imperfections of model D2 can be nat- 
urally understood as a consequence of the coarseness of the 'one- 
filament approximation' and/or uncertainties in the observations 
and the atomic data. The global upward re-calibration of the IR 
line fluxes required by this model is warranted because the gap be- 
tween the SWS and direct estimates for F(H/3) is then substantially 
narrowed. The dual-abundance models generally improve the fit to 
several lines other than the CNO ORLs, and D2 represents the most 
satisfactory compromise since the scatter in the departure ratios of 
IR CELs is much reduced. An important consequence of the above 
discussion is that it is quite plausible that the posited H-poor com- 
ponent is enriched (relative to H) not only in CNONe, but also in 
heavier elements. 

3.4.3 Ionization structure 

In Fig. 3 comparison is shown between the ionic fractions of N"*" 
and O^ in the two gas phases of model D2, while in Figs. 4 and 
5 the comparison involves the ionic fractions of the doubly ion- 
ized species of CNONe in each of the two D2 nebular compo- 
nents. The IR [O III], [Ne III] and [S III] lines are the principal 
coolants in the clumped phase and are emitted instead of the op- 
tical [O III] coolants because they are favoured by the lower pre- 
vailing Te, while the clump density (iVe = 696 cm~^) is not suffi- 
cient to severely quench most of the lines. In this respect, the model 
is in harmony with the suggestion made by T03 who, based on 
the observed deviations of O II ORL relative intensities from their 
theoretical LS-coupUng ratios, proposed that the ORL emission in 
30 Dor (and other H II regions) originated in low density regions 
and not in exceedingly dense, supposedly ionized cliunps (see also 
Ruiz et al. 2003). 

The variation of the ionic and elemental C/0, N/0 and Ne/0 
abundance ratios across the nebula was studied (Figs. 6-8). From a 
comparison with results obtained from observations, some insight 
can be gained on the best way of determining these important quan- 
tities. Here, the model C/0, N/O and Ne/O ratios (by nvmiber) are 
well approximated by the abundance ratio of their respective dou- 
bly ionized species, e.g. N/O ~ N^+/0^+ (/ N+/0+). This is true 
for the mean abundance ratios for the whole nebula, as well as for 
the ratios within the relatively H-poor phase only, and is due to the 
high excitation state of the gas. The N"'"/0"'" ratio, often used as 
an approximation to the N/O by observers who sample the optical 
[N II] and [O II] CELs only, is about two thirds the mean elemen- 
tal ratio, while in the clumps it is lower by a factor of two than 
the respective N/O of the H-poor phase. These results reproduce 
the behaviour of N/O obtained by T03 (cf. their Table 12) based on 
pure optical and pure IR forbidden-line ratios. As was also noted by 
Rubin et al. (1988), Simpson et al. (1995) and T03, in several H II 
regions IR FS line ratios yield a consistently higher N/O than opti- 
cal line ratios. According to the models shown here the IR FS lines 
of [N III] and [O iii] give a better estimate of the N/O ratio than the 
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Figure 3. Ionic fractions of .singly ionized carbon, nitrogen and oxygen in 
the H-poor and normal gas phases of model D2 vs. the thickness of the 
model filament; 'n' and 'hp' labels as in Fig. 1. 



Figure 6. The C/O elemental and ionic abundance ratios in model D2 vs. 
the thickness of model filament; curve labels marked by 'm' refer to mean 
quantities for the whole nebula, while those marked by 'Hp' refer to quan- 
tities within the H-poor phase only. 
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Figure 4. Ionic fractions of doubly ionized CNONe in the normal gas phase 

of model D2 vs. thickness. Figure 7. As in Fig. 6, for N/0. 
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Figure 5. As in Fig. 4, for the H-poor gas phase of model D2. 



Figure 8. As in Fig. 6, for Ne/O. 
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optical [N II] and [O II] CELs, both for the nebula as a whole, as 
well as for the posited, embedded H-poor clumps also (if one was 
to image the clumps in IR nitrogen and oxygen fine-structure lines 
to derive their N/0 ratios). It should be mentioned that the same 
behaviour of the C, N, and Ne ratios relative to oxygen is obtained 
for the chemically homogeneous model SI. 

3.5 On the factors from models 

Modelled mean temperatures (To,!), mean squared temperature 
fluctuation {t^) factors (Peimbert 1967) and ionic fractions Qi^^lK) 
are shown in Table 5 for the whole shell of model D2 (upper right 
panel), the H-poor component (lower left) and the normal compo- 
nent (lower right) for the elements H, He, CNONe and S; To,i and 

values were derived using, e.g. equations 3 & 4 of P03. Un- 
der the assumption for the existence of temperature fluctuations in 
a chemically homogeneous nebular medium, the elemental abun- 
dances derived from CELs overall underestimate the true abun- 
dances due to the adoption of a non-representative (overestimated) 
electron temperature for all nebular zones. Under this premise the 
heavy-element ORLs, which in most cases yield higher abundances 
than CELs, are often thought to give a more accurate measure of the 
nebular abundances since their strengths relative to H/3 are much 
less dependent on Te (and N^) than any given CEL to H/3 line 
ratio. Because of this, various workers have adopted a method- 
ology whereby abundances of elements from carbon to iron, and 
derived from CELs, are often corrected upwards. The corrections 
are estimated by adopting a representative factor which is usu- 
ally the average of those computed from the difference between 
the B aimer discontinuity and [O ill] forbidden-line temperatures 
(t%^i), and those computed from the ORL vs. CEL abundance dis- 
crepancy factor (^ADF' ^-S- from C'^+and O'^"'" ions). P03 adopted 
t^= 0.033 for 30 Dor corresponding to rms temperature fluctua- 
tions across the nebula of 18 per cent. Similar values were deduced 
from the 0^+ ADFs by Peimbert (2003) and T03. Slightly higher 
values were found from the C'^^ ADF by the same authors. Our 
modelling results however allow for an alternative interpretation of 
the observations. 

The single component, chemically homogeneous model SI 
yields t^(H+) = 0.0045, i.e. the rms Te fluctuations over the whole 
H II zone are 7 per cent (Table 5). Previous work on photoioniza- 
tion nebular models under the assumption of chemical homogeneity 
and variable gas density did not reveal factors much in excess of 
~0.01 for a wide range of ionizing colour temperatures and densi- 
ties (Kingdon and Ferland 1995; hereafter KF95). In model SI the 

factors for doubly ionized CNONe are all less than t^(H+) as 
they should be; this is a well known result (Harrington et al. 1982; 
KF95) arising from the non-coincidence of the H+ and 0^+ zones. 
KF95 warn against a direct comparison between t%^i and model- 
derived factors, such as those listed in Table 5, because the com- 
parison is valid only when the mean ionic temperatures and fac- 
tors for H"*" and O^"*" are equal. This is not exactly true even for our 
chemically homogeneous model SI where To,i(0^"'")— To,i(H+) 
= -144K and t^{0^+) = 0.69i^(H+). For the dual abundance 
model Dl, these values are respectively — 975K and 3.5t^(yi'^). 
P03 found tBai= 0.022, interpreting this result as evidence for tem- 
perature fluctuations in a chemically homogeneous nebula. This is 
not supported by model SI. It should further be noted that both 
SI and Dl models fit the Balmer discontinuity within a 3 per cent 
margin (as well as the optical [O III] lines). Only models Dl-2 
however reproduce the intensities of the CNO ORLs. For model 
D2 globally, i^(H+) = 0.022, t'^{0'^+) = 0.078 and t'^{C'^+) = 



0.14; thus the introduction of H-deficient inclusions results in the 
derivation of higher values that approach or exceed the (empir- 
ical) Even so, this verisimilitude cannot be interpreted in the 
framework of the original Te-fluctuation hypothesis since the em- 
pirical has no direct physical correspondence in the context of 
these models. This is also evident from Fig. 1 where the electron 
temperature profile across the physical extent of the two gas phases 
of model D2 is plotted. A substantial temperature gradient can only 
be discerned throughout the H-poor phase whose Tc rises, almost 
monotonically, from ~3000K at the inner edge to 6200 K at the 
outer model radius. Similarly, the empirical t adf values cannot be 
meaningfully compared with the model's values for various ions. 



3.6 Model versus empirical elemental abundances 

The abundances derived from our 4 models are presented in Ta- 
ble 6 (cols. 2-5), along with our adopted preferred estimates (col. 
6; mean of models Dl and D2) and those from previous determi- 
nations (cols. 7-10). They do not include estimates for depletion 
onto dust grains, so they pertain to the gaseous phase only. Both 
the t"^ = 0.0 and t'^ = 0.033 results of P03 are listed (the latter 
were favoured by P03). Our preferred abundances are lower overall 
than those for = 0.033, but nitrogen is 10 per cent larger; this is 
partly due to the coarse correction for unseen ions applied by P03 
(see Sect. 3.4.3). The differences are largest, about a factor of 2, for 
C and Fe. These are the elements which are mostly affected by the 
upward correction under the assumption of = 0.033. 

At this point it should be noted that abundances determined 
from detailed photoionization models and empirical methods can- 
not be immediately ranked. Models have the obvious advantage of 
self-consistency and of optimized (although possibly wrong) ion- 
ization correction factors, whereas empirical methods exactly 'fit' 
diagnostic lines by definition. The presented models exactly fit all 
the standard diagnostic Unes. While this accuracy is iUusory for ob- 
vious reasons, it does ensure that, for given atomic data, the model 
abundances arguably supersede the empirical ones. 

The ^ empirical approach not only supersedes the clas- 
sical one {t^ = 0) in that it integrates new spectral features, but 
it also takes advantage of fitting approximately the heavy element 
ORLs to claim superiority over models which do not fit them, de- 
spite the fact that there is no strong physical basis for the existence 
of large amplitude temperature fluctuations in chemically homoge- 
neous nebulae. * Here we exhibit photoionization models based on 
known standard physics, which account for the ORLs as well as for 
other optical spectroscopic diagnostics thanks to the fundamental 
assumption of chemical inhomogeneity, which we consider to be 
plausible on independent grounds (Section 4). In addition, despite 
the demanding condition of physical self-consistency, the models 
account well for a number of independent observables, while the re- 
maining moderate discrepancies are understood in terms of identi- 
fied uncertainties and the coarseness of the single filament descrip- 
tion. Therefore even though the models are probably not unique, 
their (averaged) abundances provide a basis on which empirical 
abundances can be evaluated in this particular situation. 

A comparison of the results is meaningful if differences in the 
employed atomic data and possible inaccuracies in the ionization 



^ It should be noted, however, that Stasinska & Szczerba (2001) proposed 
nebular models which resulted in enhanced values by considering small 
scale density variations in a dusty medium. 
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Table 5. Temperatures (in K), factors and mean neutral and ionic fractions, for the chemically homogeneous model SI, and the dual-abundance model D2 
(for the whole nebula and each of the two components) for H, He, CNONe, and S. The numbers in parentheses are powers of 10. 



Element 


I 


II 


III 


IV 


I 


II 


III 


IV 








SI 








D2 




H (To,i) 


9958 


9962 






10920 


9654 






it") 


1.29(-2) 


4.45(-3) 






2.09(-2) 


2.23(-2) 






(X+'/X) 


9.66(-3) 


9.90(-l) 






8.03(-3) 


9.92(-l) 






He 


10150 


9960 


9807 




9741 


9615 


9174 






1.28(-2) 


4.45(-3) 


1.65(-5) 




6.72(-2) 


2.52(-2) 


1.70(-2) 






9.35(-3) 


9.91(-1) 


1.57(-4) 




7.70(-3) 


9.92(-l) 


2.92(-4) 




C 


10860 


10610 


9904 


9592 


5421 


6585 


7821 


9089 




3.76(-3) 


5.83(-3) 


3.92(-3) 


1.03(-3) 


2.11(-1) 


2.28(-l) 


1.44(-1) 


2.84(-2) 




3.58(— 5) 


9.92(-2) 


8.66(-l) 


j.4D(— 2) 


1.06(— 4) 


1.99(-1) 


7.47(-l) 


5.31(— 2) 


N 


9433 


10790 


9876 


9600 


7162 


7129 


7901 


8789 




i.oy(— 2) 


4.68(-3) 


3.58(-3) 


1.12(-3) 


1.75(-1) 


2.07(-l) 


1.37(-1) 


5.17(-2) 




2.49(— 3) 


1.07(-1) 


8.53(-l) 


Q nit T\ 


1.60(— 3) 


1.91(-1) 


7.40(-l) 


o.oy(— 2) 


O 


9868 


10730 


9818 


9808 


9570 


7533 


8679 


8425 




1.3i(— 2) 


4.58(-3) 


3.06(-3) 




9.4D(— 2) 


1.96(-1) 


7.78(-2) 


7.37(-2) 




8.66(— 3) 


1.58(-1) 


8.33(-l) 


6.48(— 5) 


6.07(— 3) 


2.30(-l) 


7.64(-l) 


7.10(— 5) 


Ne 


9172 


10380 


9901 


9797 


5363 


5843 


8046 


8031 




Z.ZZ( — A) 


5.58(-3) 


4.02(-3) 


lAjy—j ) 


i.4U(, — I ; 


2.27(-l) 


1.30(-1) 


1 (\A ( 1 't 




1.33(— 3) 


1.28(-1) 


8.71(-1) 


D.4j( — O) 


oil/ A\ 

o.l3(— 4) 


1.44(-1) 


8.56(-l) 


1.56(— 6) 


S 


10270 


10550 


10130 


9600 


6885 


8134 


8348 


9152 




1 A'2/ 1\ 


7.64(-3) 


4.90(-3) 


1.06(— 3) 


2.38(— 1) 


1.70(-1) 


1.20(-1) 


2.4U(— 2) 




2.4d(— 0) 


3.31(-2) 


6.31(-1) 


3.21(-1) 


3.51(— 6) 


5.49(-2) 


6.72(-l) 


O CO/ 1 \ 






D2 (H-poor gas) 






D2 (Normal gas) 




H (To,i) 


5325 


4120 






10500 


9897 








3.70(-2) 


6.29(-2) 






6.71(-3) 


7.20(-3) 






(X+'/X) 


3.04(-2) 


9.70(-l) 






1.80(-2) 


9.82(-l) 






He 


5158 


4107 


3128 




10370 


9903 


9407 






4.48(-2) 


6.24(-2) 


1.51(-4) 




7.83(-3) 


7.23(-3) 


2.70(-5) 






4.79(-2) 


9.52(-l) 


2.05(-4) 




1.13(-2) 


9.88(-l) 


2.90(-4) 




C 


4699 


4672 


3847 


3421 


10740 


10550 


9850 


9465 






5.74(-2) 


5.17(-2) 


Z. Ji^ — ^) 


\jSjD\ — D) 


7.85(-3) 


6.68(-3) 


LJ\D\ — J) 




2.56(-4) 


3.80(-l) 


6.12(-1) 


7.89(-3) 


2.85(-5) 


1.26(-1) 


7.94(-l) 


1.95(-2) 


N 


5304 


4875 


3718 


3376 


10350 


10840 


9800 


9484 




3.62(-2) 


4.78(-2) 


4.10(-2) 


1.84(-2) 


8.26(-3) 


6.89(-3) 


5.96(-3) 


2.63(-3) 




1.53(-2) 


3.64(-l) 


6.01(-1) 


1.97(-2) 


2.21(-3) 


1.30(-1) 


7.78(-l) 


9.04(-2) 





5510 


4624 


3601 


3131 


10460 


10790 


9736 


9407 




2.67(-2) 


5.38(-2) 


3.45(-2) 


1.48(-4) 


6.45(-3) 


7.38(-3) 


5.35(-3) 


5.11(-5) 




2.81(-2) 


4.91(-1) 


4.81(-1) 


4.07(-5) 


1.59(-2) 


1.52(-1) 


8.32(-l) 


7.85(-5) 


Ne 


5175 


4587 


3946 


3129 


9864 


10250 


9885 


9393 




4.61(-2) 


5.93(-2) 


5.81(-2) 


2.88(-4) 


7.00(-3) 


8.11(-3) 


7.11(-3) 


3.95(-5) 




6.43(-3) 


3.17(-1) 


6.77(-l) 


8.48(-7) 


6.15(-4) 


6.26(-2) 


9.37(-l) 


1.91(-6) 


S 


4839 


4933 


4070 


3378 


10590 


10740 


10080 


9465 




5.55(-2) 


5.24(-2) 


5.95(-2) 


1.96(-2) 


6.61(-3) 


6.64(-3) 


7.71(-3) 


2.32(-3) 




1.23(-5) 


1.43(-1) 


8.05(-l) 


5.22(-2) 


3.00(-6) 


5.22(-2) 


6.16(-1) 


3.05(-l) 



correction factors adopted in the empirical approach are neutral- 
ized. This can be achieved by considering the results differentially. 
For this purpose, the = abundances obtained by P03 can be 
associated to models Sl-2 (fit to optical/UV CELs), and the = 
0.033 abundances to models Dl-2 (fit to both optical/UV CELs 
and ORLs). Excluding helium and elements not considered by 
P03, the average abundance ratios of Dl over SI and D2 over S2 
are both equal to 1.20 ± 0.10 (the aberrant N abundance of Dl 
was discarded; see discussion on [N II] in Sect. 3.4.2), while the 
same ratio for the 't^ = 0.033' over = 0' abundances is 1.48 
± 0.18. Although we have no clear explanation of why this is so, 
we note that the square root of the abundance ratio', 1.21 ± 
0.07, is identical to the 'D/S model ratio' (On close inspection 
of Table 6, however, the individual abundance enhancements 



differ significantly in the two approaches). Thus, compared to 
this particular self-consistent physical model which quantitatively 
accounts for observation, the 't^ = 0.033' abimdances are sig- 
nificant overestimates, while the mean Dl-2 model abundances 
are equal on average to the geometric mean of the t'^ = and 
non-0 empirical abundances. In other words, there is objectively no 
reason in this example to prefer one or the other set of empirical 
abundances. On the other hand, concerning physical models, it is 
noted that whereas a model which fits existing observational data 
might be wrong, one which does not almost certainly is. On these 
grounds we deem that adopting average abundances from models 
Dl-2 as representative is a justified improvement compared to 
values derived from either chemically homogeneous models or 
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Table 6. Gas-phase abundances (by number) from single and dual abundance models compared to empirical 
abundances, in units such that logAr(H) = 12. 





SI 


S2 


Dl 


D2 


This work 


POS" 


P03'^ 


T03' 


¥02^= 








Mean 


Mean 


adopted 


(2 = 0.0 


= 0.033 


= 0.0 


= 0.0 


He 


10.93 


10.93 


10.92 


10.92 


10.92 


10.94 


10.93 


10.96 


10.97 


C 


7.64 


7.66 


7.78 


7.78 


7.78 


7.80 


8.05 






N 


7.24 


7.15 


7.25 


7.24 


7.24 


7.05 


7.21 


6.89 


6.84 





8.36 


8.37 


8.45 


8.45 


8.45 


8.33 


8.54 


8.34 


8.25 


Ne 


7.64 


7.60 


7.73 


7.73 


7.73 


7.65 


7.83 


7.66 


7.84 


Mg 


6.72 


6.60 


6.74 


6.65 


6.69 










Si 


6.40 


6.45 


6.46 


6.53 


6.50 










S 


6.79 


6.83 


6.83 


6.85 


6.84 


6.84 


6.99 


6.11 


6.82 


CI 


4.77 


4.84 


4.82 


4.90 


4.87 


4.75 


4.82 


4.85 




Ar 


6.08 


6.10 


6.12 


6.16 


6.14 


6.09 


6.26 


6.15 


6.21 


Fe 


6.07 


6.04 


6.14 


6.11 


6.12 


6.25 


6.39 







« Peimbert (2003). 

Tsamis et al. (2003a). 
= Vermeij et al. (2002b). 



empirical methods. 



3.7 On the composition of the H-poor gas 

3. 7. 1 Helium and the primordial He abundance 

The question of whether He is enriched in the clumped gas along 
with CNONe is of interest, since helium is expected to be present in 
e.g. nucleosynthetic products expelled from evolved, massive WR- 
type stars (e.g. Rosa & Mathis 1987; Kobulnicky et al. 1997). The 
larger the helium abundance is in the H-deficient gas, the smaller it 
is in the 'normal' nebular component. A recent value for the primor- 
dial He/H is about 0.0786 by number (Luridiana et al. 2003). From 
model trials, the condition that He/H > 0.08 in the normal compo- 
nent translates into He/H < 0.12 for the posited H-poor clumps in 
30 Dor. Thus helium is not much enhanced relative to hydrogen in 
the clumps, and the excess heavier elements there should originate 
in extreme nucleosynthesis of the kind expected from supernova 
explosions (e.g. Woosley et al. 2002), rather than from Wolf-Rayet 
stellar winds (see Section 4). 

In both Dl-2, He/H = 0.10 , was arbitrarily adopted for the H- 
poor component (between the extremes of 0.08 and 0.12; Table 2). 
Decreasing He/H in the clumps below 0. 10, an increase of the over- 
all He/H abimdance follows imtil the latter becomes approximately 
equal to the one found from chemically homogeneous models (or 
standard empirical analyses). If, reasonably, it is assumed that He/H 
in the H-poor gas is at least as large as in the normal gas, it naturally 
follows that the He/H ratio derived from standard analyses of H ii 
regions, under the assumption of nebular chemical homogeneity, is 
only an upper limit to the actual He abundance. A similar point was 
made by Pequignot et al. (2002) in the case of PNe. This system- 
atic effect is a possible new source of uncertainty in the important 
quest for the primordial He abundance from GHIIR spectroscopy of 
blue compact galaxies (e.g. Luridiana et al. 2003; Izotov & Thuan 
2004). 

3. 7.2 Elements beyond helium 

While the major part of oxygen and other heavy elements, such 
as Ne, Mg, Si, S, CI, and Ar, is ejected into the ISM during the 



final stages of massive star evolution (M > 8M0) and the result- 
ing type n supernova explosions, the exact contribution of massive 
stars vs. low- and intermediate-mass stars in the production of C 
and N is under debate: e.g. Carigi et al. (2005) advocate Galactic 
chemical evolution models according to which only about half of 
the carbon found at present in the solar vicinity has been produced 
by massive stars. SN ejecta are strongly deficient in hydrogen and 
constitute promising sources of 'oxygen-rich clumps' such as those 
favoured in the present work (see also Section 4). Are then the 
CNO(Ne) abundances of the H-deficient clumps discussed herein 
in keeping with those expected from massive star nucleosynthesis? 
First, it must be noted that the photoionization model abundances 
for such H-poor clumps are likely to be affected by large uncer- 
tainties. Concerning nitrogen in particular, only a few weak N 11 
ORLs are available, which can result in an overestimated N/0 in 
cases of unforseen blends with even weaker lines; in addition, the 
N II lines belonging to multiplet V3 are known to be potentially 
partly excited by continuum fluorescence, not only by recombina- 
tion (see e.g. T03); this also points towards an upper limit N/O ratio 
for the H-poor gas. The weakness of the N 11 ORLs certainly im- 
plies that the N/O of the posited clumps is much less than unity, 
but it does not provide a precise value or even a secure lower limit. 
More generally, uncertainties can also arise from the fact that the 
abimdances of the climips were obtained from recombination lines 
after subtracting the contribution to their emission from the 'nor- 
mal' composition nebular gas (40-50, 50-70 and 65-70 per cent for 
C II, N II and O II ORLs respectively). Thus, while there is abso- 
lutely no doubt that the observed intensities of at least the C 11 and 
O II ORLs are larger than those obtained from chemically homoge- 
neous models,^ the exact excess, here some 55, 40 and 30 per cent 
for C II, N II and O II respectively, is not very accurately determined 
due to observational uncertainties. Another source of imcertainty in 
the resulting abundances is the coarseness of the sampling of the 
material effectively observed, given the necessarily schematic ge- 
ometry of the dual-abundance photoionization model and the fact 
that spatially-resolved observations of the posited H-poor clumps 
are not available yet; the latter goal might some day be achieved 

^ As mentioned in Section 1 this is amply documented in empirical analy- 
ses of many H II regions (Peimbert et al. 1993; Esteban et al. 2005), includ- 
ing 30 Dor (T03; P03). 
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via means of integral-field unit spectroscopy with very large tele- 
scopes. 

If we consider that our model H-poor clumps would in real- 
ity correspond to a homogeneous mixture of extremely H-deficient 
material to which 'normal' ISM gas has been partially introduced, 
the (corrected) relative abundances in the model clumps of (C/0, 
N/0) ~ (0.36, 0.09) (model D2; Table 2) can be compared to the 
(C/0, N/0) ~ (0.20, 0.03) ratios expected from typical massive star 
nucleosynthesis for the presupemova wind plus type II SN ejecta 
(cf. fig. 27 of Woosley et al. 2002). The fact that the photoioniza- 
tion models tend to indicate C/O and N/0 ratios 2-3 times larger 
than 'expected' may be a potential difficulty for either the model 
assumptions or the proposed origin of the posited clumps, and re- 
quires further analysis: the details of the mixing process of H-poor 
with H-rich gas should especially be investigated. The same con- 
clusion applies to Ne/O, although in this case only upper limits 
to the Ne ii ORLs are available and Ne/O relies on IR CELs, no- 
tably the [Ne II] 12.8-^m transition (25-30 per cent flux accuracy). 
Nonetheless, keeping in mind the aforementioned uncertainties, the 
departures between model results and massive star nucleosynthetic 
theory are not large enough to jeopardize the basic features of the 
model. The N/0 ratio is found to be definitely very small in the 
H-poor gas and therefore satisfactory at the semi-quantitative level 
currently attainable, even for cases where part of the nitrogen that 
was ejected during the presupemova evolution might not reside in 
oxygen-rich clumps. On the other hand, the determination of C/O 
in the ISM of the LMC involves the difficult calibration of the C ill] 
A1909 intensity relative to the optical, which in turn can impact the 
accuracy of the estimated contribution of the 'normal' gas to the 
emission of C ii A4267. Finally, the photoionization model analysis 
indicates that the IR lines are overall best accounted for by model 
D2 in which elements heavier than neon, such as the essentially 
better constrained sulfur and argon, are enhanced in the H-poor gas 
by similar amounts to oxygen (Sect. 3.4.2); this is expected for ma- 
terial akin to type n SN ejecta. 



4 DISCUSSION AND CONCLUSIONS 

Self-consistent photoionization models of the giant 30 Dor nebula 
in the LMC have been constructed, precisely fitting all the opti- 
cal spectroscopic plasma diagnostics available for a bright filament 
observed from UV to radio wavelengths. A dual-abimdance model 
(D2) is favoured which incorporates small scale chemical inhomo- 
geneities in the form of relatively H-deficient inclusions embed- 
ded in gas of typical LMC composition, and subjected to a CoStar 
model atmosphere ionizing spectrum. The model filament, whose 
gas pressure increases outward by a factor of 2, is essentially ra- 
diation bounded and its geometrical thickness and distance from 
the central source are in good agreement with what is suggested 
from H/3 images of the nebula. The relatively cool and dense in- 
clusions are uniformly distributed throughout the filament and are 
in local pressure equilibrium with their surroundings, representing 
2 per cent of the nebular mass, while emitting about 8 per cent 
of the total H,3 flux. This model, designed to closely fit the CNO 
ORLs (in addition to the UV/optical CELs), leads to significantly 
improved predicted intensities for a number of other optical and IR 
lines, which were not used in the convergence process. The helium 
abundance relative to hydrogen is about normal in the inclusions, 
while all other elements for which enough information is available 
(C, N, O, Ne, S, and Ar) have abundances that are enhanced by 
typically one order of magnitude compared to their respective val- 



ues in the ambient H II region. To our knowledge, this is the first 
time that an H II region model is exhibited which, incorporating 
only standard, dust-free physics for photoionized nebulae, provides 
a satisfactory solution to the extensively researched issue of dis- 
cordant ORL vs. CEL heavy-element abundances in this class of 
objects. 

The abvmdances of heavy elements in the chemically inhomo- 

geneous 30 Dor models were found to be lower than those advo- 
cated by authors who consistently apply upward corrections to the 
forbidden-line abundances of H II regions by postulating temper- 
ature fluctuations in a chemically homogeneous medium. For this 
particular nebula, the presence of H-deficient inclusions warrants 
an overall upward revision of the heavy element abundances by at 
most 0.1 dex, compared to empirical methods or models that do not 
take into account the heavy element ORLs, and indicates that the 
latter techniques provide only an upper limit to the helium abun- 
dance in H II regions. Inasmuch as, in the current paradigm, the re- 
combination lines do not provide an accurate estimate of the overall 
heavy-element to H abundance ratios in a nebula, our findings cast 
new light on the use of CNONe ORLs as tracers of abundance gra- 
dients across galactic disks such as those exhibited by Esteban et 
al. (2005). It could prove worthwhile to investigate these gradients 
further in the context of dual-abundance H II region models. 

Notwithstanding the positive aspects of the presented models, 
we must ask ourselves: is there any basis to justify the existence 
of H-deficient inclusions in 30 Dor or other H II regions and what 
could their origins be? It is often argued that since the ISM appears 
to be chemically homogeneous (e.g. Moos et al. 2002; Andre et 
al. 2003), extensive mixing of the various chemical species should 
somehow occur. However, observations establish homogeneity by 
comparing different sightlines over typical ISM depths of ~ I kpc, 
many orders of magnitude larger than the inclusions considered 
here. In fact, the 0.5 kpc local bubble appears to be oxygen defi- 
cient relative to the average few kpc neighbourhood (Andre et al. 
2003). Our interpretation of 30 Dor as a nebula which appears to be 
homogeneous on a large scale, but is chemically inhomogeneous at 
the smallest sub-parsec sized scale, resonates with the findings of 
de Avillez & Mac Low (2002) who studied the mixing of chemi- 
cal species in the ISM in the context of supernova-driven models 
and concluded that: 'even if the major fraction of the ISM gas is 
completely mixed, lower level inhomogeneities will still be seen 
as abundances are measured more carefully'. We take it that the 
recent wealth of abundance determinations in nebulae involving 
optical recombination Unes, and extensively quoted in this work, 
have indeed increased the level of accuracy with which such studies 
can be done. Those authors established that the smallest ISM mix- 
ing timescale, corresponding to the largest assumed SN explosion 
rate, is about 120Myr, with weak abundance inhomogeneities tak- 
ing longer times to be erased than strong ones due to an exponen- 
tially declining mixing rate for weaker inhomogeneities. As they 
noted however, this representative timescale is longer than the time 
interval between type II SNe and other polluting factors, such as 
e.g. H-deficient winds from evolved massive stars or star formation 
debris, which in the meantime introduce fresh inhomogeneities into 
the ISM; it is also longer than the typical lifetime of H II regions 
hke 30 Dor and their embedded stellar OB associations. Therefore, 
according to the study of de Avillez & Mac Low, the ISM may 
never reach homogeneity under such conditions. 

In an attractive scenario, Tenorio-Tagle (1996; hereafter 
TT96) proposes that, after a long journey in a hot coronal phase 
(T > 10® K), metal-rich material ejected by supemovae (SNe) 
erupting within OB associations will eventually fall back on to the 
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galactic disk as a rain of cool, dense 'droplets' of final typical size 
of ~10^^ cm and densities of ~10^ cm~^. These will be incorpo- 
rated quite evenly into the atomic H and molecular H2 content of 
the disk, without significant mixing, until a next generation of mas- 
sive stars photoionizes the gas and triggers complete mixing. Thus 
according to TT96, chemical homogeneity in a galactic disk is not 
instantaneously achieved soon after the end of a starburst event, 
but results as a consequence of the birth of the next massive star 
generation. It is tempting to identify the remains of the droplets 
singled out by TT96 with the H-deficient clumps introduced in the 
present photoionization models, especially as their size, a few per- 
cent of the size of a complete photoionized filament, and expected 
densities, ideally fits our needs. According to this depiction, one 
can envisage a galaxy's ISM abounding with small-scale chemical 
inhomogeneities which are only observed emitting metallic recom- 
bination lines, such as C 11 A4267, upon their photoionization in 
the vicinity of a young stellar cluster. In this context, we note that 
whereas the R136 cluster within 30 Dor is of young age (1-5 Myr; 
Massey & Hunter 1998) and only few SNRs have been unambigu- 
ously identified in its wider vicinity (Lazendic et al 2003; Chu et 
al. 2004), the older association Hodge 301 located 3' to the north- 
west, may have hosted ~40 SNe since its formation 20-25 M5T ago 
according to estimates by Grebel & Chu (2000). 

As emphasized by Scalo & Elmegreen (2004), a general prob- 
lem with timescale arguments pertaining to ISM mixing, is that 
they apply to turbulent transport (e.g. Roy & Kunth 1995), and not 
to chemical homogenization. For genuine mixing to be effectively 
achieved, viscosity and molecular diffusion are ultimately required, 
but these mechanisms act over very small scale lengths and/or very 
long times. Molecular diffusion becomes effective in the presence 
of strong gradients of concentration, e.g. at the interface of two 
media with very different composition. Nonetheless, the different 
thermal behaviour of such media could also delay the mixing. This 
is a very complex problem for which a generic solution has yet 
to be found. Along the same lines, the detailed simulations of de 
Avillez & Mac Low (2002) point to a very long timescale, over 
10* yr, for the final decay of small ISM inhomogeneities. Future 
work may involve a coupled hydrodynamics/photoionization com- 
putational approach to this problem (e.g. Lim & Mellema 2003). 

In the last step of the scenario developed by TT96, the dis- 
ruption of the H-poor droplets by photoionization proceeds hydro- 
dynamically aided by the presence of shocks, and the gas, now at 
10* K, acquires densities similar to those of the ambient H 11 re- 
gion; diffusion of the disrupted droplets is then suddenly acceler- 
ated. TT96 estimated a final diffusion time of the order of 10* yr in 
the lowest density H II region gas (at the small spatial scale of the 
metal-rich droplets), and a time significantly shorter than 10® yr in 
denser filaments. At first view, this timescale for complete chemical 
homogenization in the newly born H II region seems uncomfort- 
ably short to justify identification of the droplets with the H-poor 
clumps posited in the present photoionization models. Importantly 
however, we note that, unlike the supposition made by TT96, even 
in the event of its re-ionization, the H-poor gas can maintain a very 
low temperature due to fine-structure line cooling and remain rela- 
tively dense and isolated from its surroundings for some time. Thus, 
considering that the thermal properties of the H-poor inclusions 
will tend to significantly delay diffusion, it is quite possible that 
chemical inhomogeneities will survive for a significant fraction of 
the lifetime of a GHIIR. Under these conditions, identification of 
our model H-deficient clumps with the relics of T96's metal-rich 
droplets appears to be a reasonable guess. Additional work will be 
needed to develop observational tests of the model and to further 



constrain the abundances in the oxygen-rich inclusions. Whether a 
significant amount of nitrogen is present or not in the H-poor phase 
is critical with regard to the latter's origin in material lost e.g. dur- 
ing the presupemova evolution of massive stars or in a type II su- 
pernova event (Section 3.7.2). Confirming with improved accuracy 
that elements beyond neon are overabundant in the clumps in the 
same way as oxygen (Section 3.4.2), is equally important. 

In conclusion, the photoionization model description pro- 
posed here, which primarily aimed at quantitatively solving a long- 
standing spectroscopic anomaly using known, conventional H II 
region physics, also turns out to be in qualitative agreement with 
independent theoretical considerations regarding the chemical evo- 
lution of the ISM. Detailed studies of other young, relatively metal- 
poor GHIIRs, such as 30 Dor, may be a promising way to unravel 
small-scale composition inhomogeneities and provide additional 
observational support to the interesting thesis advocated by TT96. 
Efforts should be made for the detection of ORLs such as C 11 
A4267 and O II A4650 (and others of nitrogen and neon ions) in 
H II regions of all kinds, in order to build abundance determinations 
on firm grounds. Whether chemical inhomogeneity will prove to be 
the rule or the exception is an important question in its own right. 

We note that the survival of hydrogen-deficient inclusions in 
the ISM could be of paramount importance for the process of star 
formation since these inclusions, cooler and denser than their sur- 
roimdings, could act as seeds to trigger protostellar collapse. 
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